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ABSTRACT 

The c2d Spitzer Legacy project obtained images and photometry with both IRAC and MIPS instru- 
ments for five large, nearby molecular clouds. Three of the clouds were also mapped in dust continuum 
emission at 1.1 mm, and optical spectroscopy has been obtained for some clouds. This paper combines 
information drawn from studies of individual clouds into a combined and updated statistical analysis of 
star formation rates and efficiencies, numbers and lifetimes for SED classes, and clustering properties. 
Current star formation efficiencies range from 3% to 6%; if star formation continues at current rates for 
10 Myr, efficiencies could reach 15% to 30%. Star formation rates and rates per unit area vary from 
cloud to cloud; taken together, the five clouds are producing about 260 M Q of stars per Myr. The star 
formation surface density is more than an order of magnitude larger than would be predicted from the 
Kennicutt relation used in extragalactic studies, reflecting the fact that those relations apply to larger 
scales, where more diffuse matter is included in the gas surface density. Measured against the dense gas 
probed by the maps of dust continuum emission, the efficiencies are much higher, with stellar masses 
similar to masses of dense gas, and the current stock of dense cores would be exhausted in 1.8 Myr on 
average. Nonetheless, star formation is still slow compared to that expected in a free fall time, even in 
the dense cores. The derived lifetime for the Class I phase is 0.54 Myr, considerably longer than some 
estimates. Similarly, the lifetime for the Class SED class, 0.16 Myr, with the notable exception of 
the Ophiuchus cloud, is longer than early estimates. If photometry is corrected for estimated extinc- 
tion before calculating class indicators, the lifetimes drop to 0.44 Myr for Class I and to 0.10 for Class 
0. These lifetimes assume a continuous flow through the Class II phase and should be considered me- 
dian lifetimes or half-lives. Star formation is highly concentrated to regions of high extinction, and the 
youngest objects are very strongly associated with dense cores. The great majority (90%) of young stars 
lie within loose clusters with at least 35 members and a stellar density of 1 M Q pc~ 3 . Accretion at the 
sound speed from an isothermal sphere over the lifetime derived for the Class I phase could build a star 
of about 0.25 Mq, given an efficiency of 0.3. Building larger mass stars by using higher mass accretion 
rates could be problematic, as our data confirm and aggravate the "luminosity problem" for protostars. 
At a given Tboi, the values for L^oi are mostly less than predicted by standard infall models and scatter 
over several orders of magnitude. These results strongly suggest that accretion is time variable, with 
prolonged periods of very low accretion. Based on a very simple model and this sample of sources, half 
the mass of a star would be accreted during only 7% of the Class I lifetime, as represented by the eight 
most luminous objects. 

Subject headings: stars: formation — infrared: stars — ISM: dust — ISM: clouds 
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1. INTRODUCTION 

Star formation in the solar neighborhood has been stud- 
ied extensively, but a complete survey for forming stars 
in nearby molecular clouds has been lacking. The IRAS 
catalogs provided an unbiased survey, but they were lim- 
ited in sensitivity, wavelength coverage, and spatial reso- 
lution. For objects in embedded phases of evolution, the 
luminosity limit w as about 0.1(d/140) 2 L with d m pc 
(My ers et al.lll987l ). The wavelength coverage from 12 to 
100 fim covered the peak of emission from typical Class I 
objects, but deeply embedded Class sources were often 
missed, requiring special processing to be dete cted in IRAS 
data after being found at longer wavelengths ([Andre et alj 
2000). In addition, shorter wavelength data were needed 
to identify less embedded objects and to learn more about 
the embedded ones. Follow-up studies at shorter wave- 
lengths often found multiple sources, which could be con- 
fused within the large beams of IRAS. Detailed studies in 
the mid-infrared from the ground were hampered by the 
atmosphere and limited to small regions already known to 
have luminous sources or clusters. 

Ground-based near-infrared surveys, enabled by the 
development of large format near- infrared arrays, are a 
powerful tool for finding low luminosity YSOs in nearby 
clouds with relatively low extinction. High spatial resolu- 
tion and sensitivity in these surveys d ecrease the level of 
confusion found in older surveys (e.g., IStrom et al]|1989l 
lEiroa et al.lll992t) . Indeed, the 2MASS survey at J, H 
and K has been used extensively to characterize the young 
star population and d istribution in a number of clouds 
(e.g.. iCarpenter ||2000| ). Many of these studies have fo- 
cused on young clusters in m ore massive star-for ming re- 
gions, however (see review bv lLada fc Ladall2003D . 

The Infrared Space Observatory (ISO) improved signif- 
icantly on IRAS in terms of sensitivity and spatial res- 
olution at the mid-infrared wavelengths. The ISOCAM 
instrument mapped a few square degrees of nearby star- 
forming clouds in two bands at 6.7 and 14.3 /im (e.; 
IBontemps et al.ll200lL iPersi et al.ll2003l . iKaas et alj|200l 
The mid-infrared observations revealed the more embed- 
ded population compared with near-infrared surveys. The 
number of Class II sources with Lh n i > -03 L Q in Ophi- 
uchus was doubled (|Bontemps et al.ll2001h . These studies 
were focused on the densest parts of the clouds, but they 
clearly indicated the importance of deep mid-infrared ob- 
servations. 

The Spitzer mission, combined with ground-based near- 
infrared and submillimeter data, can remedy most of the 
deficiencies of previous studies. Spitzer can provide a 
survey that is deep, wide, and relatively unbiased, and 
that covers the entire wavelength range needed to charac- 
terize the different evolutionary stages, from deeply em- 
bedded YSOs to young stars which have lost most of 
their disks. Here we summarize the results of the Spitzer 
legacy project, "From Molecular Cores to Planet-forming 
Disks", or "Cores to Disks" , further abbreviated to c2d 
(|Evans et al.l 120031 ) . One of the main approaches of the 
c2d project has been to provide a more complete, less bi- 
ased sample of star formation in nearby large clouds and 
small cores. Toward that end, we used Spitzer to map 15.5 
square degrees in five large, nearby molecular clouds and 
about 0.6 square degrees in 82 small dense cores. For this 



paper, we focus on the studies of large clouds, where the 
mapping efficiency of Spitzer enables coverage of large ar- 
eas, in contrast to previous surveys, which focused on small 
areas around IRAS sources, for exa mple. The results o n 
the small cores will be summarized bv lHuard et alj ( 



Based on all of our data and some auxiliary data, we 
construct a combined list of young stellar or substcllar ob- 
jects (hereafter YSOs). We discuss issues of contamination 
and completeness, and we calculate various quantities that 
will be used in the analysis (SJ2]). We summarize the star 
formation efficiencies and rates (Q, comparing our results 
to predictions used in extragalacti c work fEj4.ll). Thes e val- 
ues supercede preliminary values (jEvans et al.ll2u"0 8). We 
compare numbers of sources in various SED classes and 
calculate lifetimes for these classes (31]). We then reexam- 
ine the issues of source classification (*JS]), the connection 
between empirical classes and various stages of star forma- 
tion, and the estimation of lifetimes for these stages ({2}. 
We discuss the spatial distribution of sources and their 
clustering properties (®. Finally, we compare the obser- 
vational results to the predictions of various theoretical 
models (33), describe future work f UlOp , and summarize 
the main results ( tjll[) . 

1.1. The Sample 

Fiv e large clouds wer e selected for the c2d proje ct: Ser- 
pens (lEiroa et al.ll2008f). Pe rseus dBallv et al.ll2008D. Qphi- 
uchus (IWilking et al J 120081). Lupus (|Comer6nll2008f ). and 



Chamaeleon (jLuhmanl 120081) . More specifically, we tar- 
geted Lupus I, III, and IV and Chamaeleon II (hereafter 
Cha II). These were chosen to lie within about 300 pc of 
the Sun, to span a range of previously-known star forma- 
tion activity, and to complement observations of smaller 
regions in these clouds obtained by Guaranteed Time Ob- 
servers (GTOs). We have included data from the GTO 
observations with our data for a complete picture of the 
clouds. 

The clouds in this study are listed in Table [TJ with 
the adopted distances, solid angles and areas with 
both IRAC and MIPS data, a measure of turbulence, 
the total mass in the mapped area, and the cross- 
ing time. Some of these clouds were targets of ISO 
survey s of smaller regions: 0.7 sq. deg. in Ophi- 
uchus (IBontemps et alj|200ll ); 0.13 sq. deg. in Serpens 
(IKaas et al.ll2004l) : and 0.2 sq. deg. in Cha II (jPersi et alj 
2003). The results presented here are base d on detailed 
studi e s of the individual clouds: Serpens (lHarvev et al.l 
12001 lHarvev et alj 120071 lHarvev elaTl l2007bl); Perseus 



Ophiuchus dPadgett et al.l 


2008b]: |Allen et al.ll2008b1 : Lu- 


pus (Chapman et al.ll2007: 


Merin et all2008al): and Cha II 


(lYoung et al.l 20051 Porras et all 2007; Alcala et all 12008; 



papers have been updated based on a consistent analysis 
of the combined dataset. 

None of these clouds are forming very massive stars; the 
most massive star in the close vicinity of a cloud mapped 
by c2d is HD147889, a B2 (Ill/TV) star lying within our 
map of Ophiuchus. The a Sco group is farther away (pro- 
jected distance from the cloud of 4 pc), consisting of an 
Q9V star and a B2III star and at least two other B stars 
(Pigulski 1992])- The star formation patterns in Ophi- 
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uchus may have been affected by in teraction with these 
stars (|Lorenlll989l . iNutter et al.ll2006t) . For Perseus, there 
are hints of an interaction with 40 Per, a BO. 5 star in the 
Per OB Associa tion, which is 26 p c from the L1451 core, 
as dis cussed by IWalawender et all (|2004h and iKirk et all 
(2006}). Serpens contains a Herbig Ae star, VV Ser, which 
produces a large nebulosit y (green object toward the s outh 
in Figure CD) discus s ed by Pontoppidan et all ((2007a) and 
iPontoppidan et all (|2007b ). In addi tion, Serpens contai n 
3 other A- type stars and a B8 star ()01iveira et al.|[2008l) . 
The most ma ssive star in the Cha II cl oud is the F-type 
star DK Cha (jHughes fc Hartigan 19921); all th e rest have 
spectral type later than K2 ( Spezzi et al.ll200lj[) . The Lu- 
pus clouds are in the vicinity of the Scorpius-Centaurus 
association, which flanks the Lupus clouds at a distance 
of about 5 degrees, or 17 pc at a distance of 200 pc. 
The sub-groups of the association are called Upper Scor- 
ius (5-6 Myrs) and Upper Centaurus-Lupus (14 Myrs) 
de Geus et aTlll989D . The strong high -energy radiation 
and supernova remnants from those OB stars might have 
had an important role in the formation and evolution of 
the Lupus clouds (|Tachihara et al.|[2001[ ) . Specificially, the 
highly fragmented cloud structure is thought to be related 
to the effects of the nearby OB association. 

There are distance uncertainties for all five clouds, which 
are discussed in the detailed papers on those clouds. Here 
we supply our adopted uncertainties and propagate those 
into the areas and masses. Recent parallax observations of 
radio emission from y oung stars in Ophiu chus yield a dis- 
tance of 120 ± 4.5 pc (jLoinard et al.ll2008l ). By combining 
exti nction maps with parallaxes from Hipparcos and Ty- 
cho, MamaiekTl 20081) f ound a distance of 135 ± 8 pc, while 



lLombardi et al I (|2008f ) have derived a distance of 119 ± 6 
pc. These are all consistent with our assumed distance and 
uncertainty of 125 ± 25 pc. Because it is not clear that the 
whole cloud is at the same distance (jLoinard et al.l [20081 
lLombardi et al.|[2008f ) , we retain our value and uncertainty 
in the analysis. A recent astrometric measurement of wa- 
ter masers in N GC1333, in Perseu s, provides a distance 
of 235 ± 18 pc (|Hirota et a l. 2008). Distance estimates 
for IC348 are larger, so we retain our standard distance 
and un certainty of 250 ±50 for the Perseus cloud as a 
whole. lLombardi et alj (|2008l ) also derived a distance to 
the Lupus complex of 155 ± 8 pc, with some evidence that 
individual Lupus clouds were at different distances, both 
consistent with our assumed values of 150 ± 20 pc for Lu- 
pus I and IV and 200 ± 20 pc for Lupus HI. We use a 
distance to Serpens of 260 ± 10 pc (jStraizvs et al.| [l996L 
as was used in all our previous papers on this cloud, but 
we note that t here is some recen t evidence for a distance 
of 230 ± 20 pc (jEiroa et al.ll2008h . 

Seeking a common measure of the turbulence of these 
clouds, we chose to give the FWHM linewidth of the 13 CO 
J = 1 — > line, ideally averaged over a full map of the 
cloud. This was not available in all cases, so we have re- 
sorted to accepting stat ements like "a typical linewidth 
is ..." (|Hara et al.l 1999T ) for Lupus. For Perseus, Ser- 
pens, and Ophiuchus, J. E. Pineda (pers. comm. 2008) 
has provided the mean and median linewidth s for the en- 
semb le of spectra in the COMPLETE maps (jRidge et al.l 
l2006f) . and the average and difference of the median and 



mean are given in Table [T] as the value and uncertainty. 
iTothill et ail (|2008l ) have provided the same kind of infor- 
mation for the J = 2 — > 1 line of 13 CO in Lupus, based 
on maps with AST/RO. The linewidths are based on di- 
viding the integrated intensity by the peak intensity and 
correcting for the channel width of 0.91 km s , assuming 
Gaussian lines. The results are substantially larger than 
the estimates from the 13 CO J = 1 — » lines. For Lupus 
I, (Av) = 2.17 ± 0.05; for Lupus III, (Av) = 2.11 ± 0.05; 
and for Lupus IV (Av) = 1.53±0.03. It would be valuable 
to have more consistent and sophisticated measures of the 
turbulence in these clouds. 

The cloud mass (gas and dust) has been derived 
from extinction maps made from our c2d d ata, together 
with 2MASS data, toward background stars (| Evans et alj 
200 7J). We used the extin ction law with Ry = 5.5 
(|Weingartner &: Draind 120011 ) , which repr oduces reason- 
ably well the data for mole cular clouds (|Flahertv et alj 
l2007t IChapman et~aTl I2009D . to determine Ay. The 
mass was then ca lculated usin g the relation Nh / Ay = 
[IMQCextiV))- 1 (|Drainel l 2003l) and the value of C ext {V) 
from the on-line tables I for Ry = 5.5. This grain 
model results in a conversion from extinction to hy- 
drogen column density of 1.37 x 10 21 cm -2 mag -1 , in- 
stead of t he usual 1.87 x 10 21 , established for diffuse 
ISM gas dBohlin et al.l Il978l). In ou r earlier papers 
(iHarvev et alj l2007bl. lAlcala et al.1 12008L IChapman et alj 



120071 iMerin et al.ll2008al) . the value for the diffuse gas was 
used, resulting in cloud masses overestimated by a factor 
of 1.4 relative to those in Table [TJ 

The mass includes all cloud mass above an extinction 
contour with Ay = 2 mag, with no assumptions needed 
about geometry. The mass derived from our extinction 
maps refers to the same area covered in our surveys, mak- 
ing it ideal for calculations of efficiency, etc. The two ex- 
ceptions to this statement are Serpens, for which our sur- 
vey was designed to cover completely down to Ay = 6, and 
Ophiuchus, for which we tried to cover completely down 
to Ay = 3. The cov erages were based on extinction maps 
bv lCambresvl ([1999D . but some areas down to Ay — 2 were 
covered incompletely. The mass inside the Ay = 6 contour 
for Serpens is 1532 M©, 76% of the mass within Ay = 2. 
The mass inside the Ay — 3 contour for Ophiuchus is 1914 
Mq, 88% of the mass within Ay = 2. We use the value for 
Ay = 2 in what follows, but the difference will not be large 
for these clouds. Table Q] also gives the cloud crossing time 
in Myr, calculated from the area and the mean speed of 
turbulent motions, approximated by (v) — 0.68Au, with 
Av the 13 CO J = 1 -v linewidth. 

With the sensitivity of Spitzer, we can detect sources 
with luminosity as low as 10 -3 L Q at 350 pc, but contam- 
ination by other kinds of sources ultimately limits our sen- 
sitivity, as discussed bel ow. A similar study of the Taurus 
cloud has been done bv lPadgett et al.l (|2008aD . and other 
clouds in th e Gould Belt are be ing mapped with the same 
techniques (|Allen et al.l 12008a *) . When those studies are 
completed, the analysis in this paper can be extended to 
cover most star formation in large clouds within about 300 
pc of the Sun. 

2. OBSERVATIONS 



Available at http://www.astro.princeton.edu^draine/dust/dust.html 
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The observations here were all obtained by the c2d 
project or by GTO observations that we have included 
in our data. They have been described in the publications 
given in Nl.ll Here, we give a brief summary of the data 
used for this paper; a more complete description can be 
found in lEvans et alj (|2007ft . 

The Spitzer instruments were used to obtain data from 
3.6 to 160 /tm, but the 160 /mi data are limited by sat- 
uration, incomplete coverage, and the large beam. Pho- 
tometry at 160 /tm is not available as a standard product 
in the c2d catalogs. We provide and use some limited 
photometry at 160 /mi here. Data from the four IRAC 
bands (3.6, 4.5, 5.6, and 8. /tin) and data f rom the 2 /tm 
All Sky Survey (2MASS) (|Cutri et al]l200!h were merged 
using a 2" matching radius. These sources were merged 
with photometry from MIPS at 24 /im (4" radius). The 
70 /im data were merged with an 8" radius, but human 
judgment was used when multiple candidates for merging 
were available at the shorter wavelengths. The merging 
of IRAC and MIPS data could be done only for the area 
covered by both IRAC and MIPS observations, and larger 
areas were usually covered by MIPS observations. Since 
we use data from both instruments to separate YSOs from 
other sources, we focus here on the areas with both IRAC 
and MIPS data. 

For one clou d (Cha II), optica l photometry has also 
been published (|Spezzi et alj 12007). Also for Cha II and 
Serpens, we have reasonabl y co mplete spectral type infor- 
mation (|Spezzi et alj 120081 and lOliveira et all I2008D . For 
three clouds (Perseus, Ophiuchus, and Serpens), com- 
plete maps of dust continuum emission at 1.1 mm were 
made using Bolocam on the Caltech Submillimeter Ob- 
servatory. These have been pub l ished (|Enoch et alj|200d 
lYoung et all 120061 lEnoch et al] 120071) and the Bolocam 
data have been combined with Spitzer data (|Enoch et al.l 
2008a, E noch et alj | 200 8b). We will draw on those results 
here. We al so have less sensitiv e maps of 1.3 mm emission 
for Cha II (lYoung et all I2005D a nd maps of selected re- 
gions at 350 /im ( Wu et aLll2007f ). Maps of emission from 
the J = 1 — » transitions of CO and 13 CO were obtained 
by the C OMPLETE proje ct for Perseus, Ophiuchus, and 
Serpens (jRidge et al.ll2006D . Maps of the J = 2 -> 1 line of 
^ CO were obtained for Lupus and Cha II bv lTothill et alj 

41x1). 

3. RESULTS 

The fundamental results for the large clouds are sum- 
marized in Figures [T] and O and Tables [1] to [7] 

Figure Q] presents three color images of all the clouds 
on the same physical scale, with a 3 pc scale bar. Fig- 
ure [2] shows the extinction in grayscale, with all clouds on 
the same physical scale and extinction grayscale. All the 
identified YSOs are shown, color-coded by their SED class 
( t)5.ip . The total area covered is dominated by Perseus, 
Ophiuchus, and Serpens. 

Tables IA8I to IA12I are described in the Appendix. They 
are electronic tables providing a full list of YSOs in the 
five clouds, based on an updated, uniform analysis of the 
entire sample, using the final data release from the c2d 
projeci[3- The information in the tables in this paper is 
a condensed version of that found in the full c2d catalogs, 

16 Available at http://ssc.spitzer.caltech.edu/legacy/c2dhistory.html 



along with supplementary information. Briefly, the c2d 
catalogs provide flux densities, uncertainties, and various 
flags for each wavelength from 3.6 to 70 /im, similar infor- 
mation for the 2MASS sources that match our sources, a 
source type, and a spectral index. Since the source type 
plays an important role in discussing contamination, we 
briefly summarize the nomenclature here. 

Only sources with detections in at least three bands 
could be classified at all. If an object was consistent with 
a (possibly reddened) stellar photosphere, it was labeled 
"star". Candidates for YSOs required detections in all 
four IRAC bands and MIPS-1; those meeting stringent 
criteria (see 13. ip are labeled YSOc, with the "c" empha- 
sizing that they are only candidates. They may have 
one or more appended suffixes, such as "red" , "PAH_em" , 
or "star+dust" . Sources labeled "red" have flux densities 
at 24 /tm that are at least three times the flux density 
from the nearest available IRAC band. Sources labeled 
"PAH_em" have colors indicative of a peak in the 8 /tm 
band. The "star+dust" designation indicates that the 
SED is consistent with that of a stellar photosphere for 
wavelengths shorter than a particular band but an excess 
in that band of at least 3 er. This band was appended to 
make a full name such as "YSOc_star+dust(MPl)" where 
MP1 implies the first MIPS band, at 24 /tm. As discussed 
below, some sources lacked photometry at enough wave- 
lengths to be classified as YSOc, but were later added as 
YSOs based on other information. These may have source 
types in Table [A12l of "rising" , in which the available flux 
densities rise significantly to longer wavelengths (15 cases), 
or "redl" , in which the source is detected only in IRAC-4 
or MIPS-1 (3 cases). 

In this paper, we refine and standardize the analysis of 
the combined data for the five large clouds, discussing is- 
sues of contamination by background sources ( ij3.ip and 
completeness f ^3.2|) . After describing the standard classi- 
fication tools provided in the catalogs, we describe analysis 
that goes beyond that provided in the catalog. We discuss 
extinction corrections fi J3.3[) that can be be made to the 
flux densities and the calculation of bolomctric luminos- 
ity (Lboi) and bolometric temperature (Tboi)- We present 
in the Appendix the list of YSOs, including flux densities 
from 2MASS (1.25 - 2.17 /mi) and Spitzer (3.6 - 70 /tm). 
Additional flux densities from observations with other tele- 
scopes, at wavelengths ranging from 0.36 /tm to 1.3 mm, 
are presented in Tables lATHlAlJIlATOl and En] This list of 
YSOs is improved over the list of YSOc supplied with our 
delivery because we have removed some suspect sources, 
added known sources, added data at other wavelengths, 
and calculated additional quantities, provided in Table 
IA12I We describe this process in the following subsections. 

3.1. Contamination by Other Sources 

The main challenge is to separate YSOs from contami- 
nants. By far, most contaminants are stars, mostly back- 
ground, without infrared excess. These have colors very 
close t o zero in a colo r-color diagram using the IRAC 
bands (|Allen et al.ll2004T ) and are easily removed. The pri- 
mary remaining contaminants are then background galax- 
ies with active star formation; these may have colors very 
similar to those of embedded young objects, so we must 

with a thorough description (Evans ct al. 2007) 
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use magnitude information as well as colors. 

The automated criteria described by I Harvey et al.l 
(2007b) were used for all five clouds to produce catalogs of 
YSO candidates (YSOc), which are part of our final data 
delivery. Compared to earlier catalogs, ad vances in dis- 
tinguis hing YSOs from background galaxies (jHarvev et alj 
20073) have provided much cleaner samples. These crite- 
ria assign to each source with sufficient information an 
unnormalized probability (P(Galc)) that the source is a 
galaxy. The calcu lation of the p r obabili ty is complex and 
fully d escribed in I Harvey et al.l (|2007bf ) and lEvans et al.l 
(2007), so we summarize it here. The probability depends 
on the location of the source in three color-magnitude di- 
agrams, whether the source is extended, and whether the 
source is detected above a threshold at 70 /jm. The as- 
signment of probabilities is entirely empirical, based on 
excluding objects in a suitable sample of extragalactic ob- 
servations. 

By comparing to a set of sources fro m the SWIRE surve y 
of the ELAIS Nl extragalactic field (jSurace et al.ll2004h . 
suitably processed to simulate the c2d sensitivity and the 
extinction distribution of each cloud, we define a thresh- 
old in P(Galc) below which a source becomes a YSOc. 
For the catalogs of the five large clouds, the aggregate 
plot is shown in Figure [3] As was true for the individual 
clouds, most YSOc separate nicely from galaxy candidates 
(Gale), but there are some sources in a lower probability 
tail from the big peak in the Gale sector that are am- 
biguous. Thus, confusion with star-forming galaxies can 
be a problem for both contamination and completeness. 
Sources with \ogP(Galc) < —1.47 were assigned the YSOc 
moniker in the c2d catalogs, b ased on the analysis of Ser- 
pens by Ha rvey et alj (|2007b| ). From the surface density 
of objects from the degraded SWIRE data that would be 
misclassified as YSOc for each cloud, multiplied by the sur- 
face area of the cloud, we estimate that there could be as 
many as 51 contaminating galaxies in our total 15.5 square 
degrees. That would still be a small fraction of the total 
of 1086 YSOc (see below). We next describe attempts to 
minimize the number of contaminants. 

Certain categories of sources were checked by eye for all 
the clouds in an attempt to eliminate residual contami- 
nants. Typically these were sources that had been "band- 
filled" at 24 fim or sources with peculiar SEDs. Sources 
that were well detected at some wavelengths received flux 
estimates at other wavelengths, based on the known po- 
sition, in the process of band-filling. Some of these flux 
estimates turned out to be emission from the wings of 
nearby bright sources and other artifacts. In addition, 
objects with the "PAH_em" suffix received extra scrutiny 
because this feature is common in star-forming galaxies. 
Objects that had galactic morphologies were removed from 
the final YSO list; however, many star-forming galaxies are 
point-like to Spitzer. A total of 91 YSOc were eliminated 
by this process. We do not know how many point-like 
galaxies were removed by this process, but the numbers 
are consistent with removing most of them. 

While stars without infrared excess are readily removed 
from our sample via color criteria, post-main-sequence 
stars with circumstel lar shells can masquerade as YSOs. 
lHarvev et all (|2007bD were able to identify and remove 4 
of these behind the Serpens cloud. Follow-up optical spec- 



troscopy toward Serpens (jOliveira et al.|[2008f ) have pro- 
vided spectral types for 78 objects, 58 of which had been 
identified as YSOs. Based on their analysis, we rejected 
11 of those 58 sources ( about 20%) as back ground giants 
with infrared excesses (lOliveira et all l2008h . All 11 are 
classified as either Class II or Class III sources, with 9 out 
of 11 (about 80%) classified as Class III. We have removed 
these 11 sources from the final sample. We so far lack the 
data required to identify and remove background giants in 
Lupus, Perseus, and Ophiuchus. However, Serpens should 
be the worst case because it lies at low galactic latitude 
and longitude. Lupus and Ophiuchus may also be affected, 
though less so than Serpens. Additionally, a spectroscopic 
study of Cha II found that 96% of the YSOs identified i n 
our catalogs were true cloud members (jSpezzi et al.ll2008h . 
These points, combined with the fact that most of the 
background giants (80%) are Class III sources, to which 
we arc incomplete anyway f ^3.2p , lead us to conclude that 
contamination by background giants will not significantly 
affect our main results. 

Adding these 11 background giants in Serpens to the 
91 sources removed as described above gives a total of 
102 sources removed from the initial sample of 1086 YSOc 
(Table [2]). 

3.2. Completeness 

Since the vast majority of sources in the catalogs are 
background stars, we require an infrared excess for a source 
to be considered a candidate YSO. Consequently, our sam- 
ple is missing Pre-Main-Sequence Stars (PMS) that no 
longer have infrared excess but may have Ha emission, 
X-ray emission, etc. Complete surveys for such objects 
are needed to complete the sample of PMS objects, so we 
concentrate on YSOs, defined here to have an infrared 
excess. More complete samples of PMS objects exist for 
some clouds that have the needed comp l ement ary data, 
but not for all. For example. lAlcala et~a l. (2008) found 51 
certain and 62 likely PMS stars in Cha II, more than dou- 
ble the number (26) of YSOs. The ratio of PMS to YSOs 
in our catalog is 2 for Class II sources and 4.8 for Class III 
sources. Cha II may be a worst case because many of the 
PMS stars are not in the region studied by c2d. Restrict- 
ing attention to that region, only 4 PMS stars are added to 
the YSOs found by c2d. A similar study in Lupus would 
add 19 PMS stars w ithout infrared excess to the 94 YSOs 
(jMerin et al.ll2008al ). 

While objects of very low luminosity may still be lost 
among the galaxy background, we estimate that our YSO 
sample is 90% complete down to a luminosity integrated 
from 1 t o 30 iim of 0.05 L w and 50% complete down to 
0.01 L Q (jHarvev et alj|2007bh for our most distant cloud 
(Serpens at 260 pc, which also has the highest density of 
background stars). For comparison, the survey of Serpens 
with ISO reached a limit of 0.08 L Q , but had only tw o 
bands and covered only 0.13 sq. deg. (|Kaas et al.ll2004l ). 
They found 61 YSOs compared to our 227, mostly because 
of their smaller areal coverage. The 1-30 /im wavelength 
range covers the bulk of emission from YSOs with infrared 
excesses arising from circumstellar disks; thus these com- 
pleteness limits are a good proxy for our completeness to 
such objects. A 2 Myr old object at the mass boundary 
between stars and brown dwarfs has a luminosity of about 
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0.01 L . Indeed some objects believed to be brown dwarfs 
with disks based on more complete analysis do not make 
our YSOc list, either because their fluxes are not of high 
enough quality or because they have \ogP{Galc) > —1.47. 
For example, 2 of the 4 added PMS objects i n Cha II are 
very low mass brown dwarfs with excesses (Eh ers et al.l 
I2007T) . 

For younger YSOs still embedded within their dense 
cores, this 1—30 /xm wavelength range is less appropriate 
since the bulk of the emission is reprocessed by the sur- 
rounding envelope to the far-infrared. A separate search 
for embedded objects with luminosities less than 1 L 
found that we reach a similar completeness limit for the 
internal luminosity of these objects {Lj nt of 4 x 10~ 3 Lp 
at 140 pc, or ~ 0.014 L at 260 pc) l|Dunham et al.ll200aT . 
where this limit is set by the sensitivity of the 70 /mi MIPS- 
2 observations. The internal luminosity measures the con- 
tribution from the embedded source, after correction for 
the effects of heating by the interstellar radiation field. 

We require good quality detections in all 4 IRAC bands 
and MIPS-1 for a source to be classified as a YSOc. This 
requirement can result in our missing two types of sources: 
deeply embedded sources that were not detected in all 
bands and strong sources that saturated the detectors. 
Note that the latter can include both embedded objects 
and more evolved YSOs no longer embedded within their 
dense cores. A total of 40 sources that are clearly YSOs 
but which did not make the automatically-generated YSOc 
list for one of the two above reasons were added by hand 
to our final sample, bringing the final sample size to 1024. 

Of these 40 sources, 36 were added by co mparison to the 
search e s for embedded objec t s presented by Dunham et alj 
(l2008fl.lEnoch et all (l2008al) . Uorgensen et all (|2007[ ). and 
Uprgensen et al.1 (|2008T ). Two were added by being well- 
known YSOs that were saturated in one or more of the 
Spitzer bands but not included in the samples of embed- 
ded objects compiled by the above authors. For all satu- 
rated sources, data from other telescopes were substituted 
for the saturated Spitzer data. Three of the added sources 
actually have source types of Gale, one each in Lupus, 
Perseus, and Serpen s. The source in Lupus was added by 
iMerin et al.l (|2008al ) because it lies essentially on the bor- 
der between YSOc and G ale, and it is cand idate to be a 
brown dwarf with a disk (| Alters et alj [20061. The source 
in Perseus was added because it is contained in the sampl e 
of embedded objects compiled by Uorgensen et ail (|2007t ). 
The Serpens source is associated with an outflow. 

Turning the luminosity completeness limit into a limit 
on stellar mass requires further analysis. For the early, 
embedded stages, the luminosity depends on the product 
of stellar mass and mass accretion rate. As discussed in 
^9.21 accretion at the mean rate expected in a Shu-type 
model would predict a luminosity of 1.6 L for a central 
object of 0.08 M . Our luminosity limit of 0.014 L would 
then translate into a mass limit of 7 x 10~ 4 M . However, 
there is strong evidence for highly variable accretion rates 
( ^9.2p . so this limit is highly suspect. For the later stages 
when accretion luminosity is negligible, spectral types and 
accurate extinction corrections are necessary to determine 
masses. For the most distant cloud (Serpens), a main- 
sequence luminosity of 0.05 L implies a mass of 0.08 
M , while 0.01 L correponds to 0.04 M (|Chabrier et alj 



2000). At young ages, the mass limi ts would be lower. In 
Serpens, the masses determined by lOliveira et al.l |2008) 
from spectral types range from 0.2 to 3.0 M . In Lupus, 
masse s of YSOs are complete down to 0.1 M ( Merm et alj 
l2008ah . and in ChalL. masses extend well below the stellar 
limit to 0.015 M (jSpezzi et al.ll2008[) . We do clearly miss 
some substellar objects with disks, as discussed above. Un- 
til spectroscopy is available for a larger fraction of the 
young objects, we can only estimate our mass complete- 
ness to be near the stellar/substellar boundary. 

3.3. Classification 

For each source with sufficient data, the c2d catalogs 
provide a least squares fit to all photometry between 2 /im 
and 24 /im in the c2d catalog to determine the spectral 
index a: 

_ dlog(XS(X)) 
a dlog{\) 1 ' 

where A is the wavelength and S{\) is the flux density 
at that wavelength. These were used to class i fy obj ects 
into the four classes defined by I Greene et al.l (I1994TI . as 
described in detail in ij5.ll Class I, Flat, Class II, and 
Class III. 

An alternative classification scheme, used especially for 
more embedded objects, employs the bolometric temper- 
ature (Tboi). We will also need the bolometric luminos- 
ity (Lf,oi) for later analysis. These were computed, using 
both Spitzer data and auxiliary data at both shorter and 
longer wavel engths, by calcula ting the first two moments 
of the SED (|Chen et al.lll995[) using integ r ation methods 
described in more detail bv lPunham et all ([2008). Uncer- 
tainties in Lboi and Tf, i are dominated for early phases 
by incomplete sampling of the SED at wavelengths from 
70 to 1000 /mi; they were estimated t o be 20 to 60% 
([Enoch et al.l l2008a. Du nham et al.l l2008). The poor sam- 
pling at far-infrared wavelengths arises in part because of 
the large beam of Spitzer, causing confusion in crowded 
regions, and from the lack of data at 350 /im. A recent 
analysis of the Serpens B cluster using the fine-scale 70 iim 
mode and adding 350 /im data (|Harvev & Dunhaml l2009) 
was able to separate sources that were confused in the c2d 
data. They find that the values of L^oi and Tf, ; change by 
about the amount estimated above. Nonetheless, only one 
source out of 18 would change classification (from Class I 
to Class 0). 

For some of the analysis, it is desirable to correct the flux 
densities for extinction. Based only on the c2d data, such 
corrections would be highly uncertain. Optical follow- 
up studies have provided spectral types for the Class II 
and Class III population s in Chamaeleon II (|Spezzi et al.l 
l2007t ISpezzi et al1l2008D Lupus (F. Comeron et al., in 
preparation; Merm et all l2008af) . and part of Serpens 
(Oli veira et al.ll2008f ). For Ophiuchus and Perseus, along 
with the remainder of Serpens, we assume a spectral type 
of K7 for the Class II and Class III sources. We then used 
the spectral type (known or assumed), the near-infrared 
photometry (J, H, and Kin order of preference) , and the 
Wcingar tner fc Draind (|200l[ ) extinction law for Ry = 5.5 
to correct the photometry for extinction. The extinction 
law was chosen to match that used for determining our 
extinction maps and cloud masses ( tjl.l|) . For Class I and 
Flat objects, we assume the mean extinction to all the 
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Class II objects in the same cloud. The idea is that this 
correction removes foreground extinction, but not local ex- 
tinction from the surrounding envelope, which will be rera- 
diated in the far-infrared. The mean extinctions to Class 
II sources for the five clouds are Ay = 3.95 mag in Cha II, 
Ay — 2.91 mag in Lupus, Ay = 5.92 mag in Perseus, 
Ay = 9.57 mag in Serpens, and Ay = 9.76 mag in Ophi- 
uchus. Extinction values calculated for Class II and III 
sources without spectral type information should be re- 
garded as highly uncertain. Extinction values for Class I 
and Flat sources should be regarded as averages for each 
cloud only and not as actual extinctions towards each ob- 
ject. 

After these extinction corrections were made, the val- 
ues of spectral index, X^, and Lt, i were recomputed and 
these values are distinguished from the observed values 
by primes (a', T b ' oi , and L' bol ). Both sets of values are 
given, along with the value of Ay used for the extinction 
correction, in Table [A"T2l When wavelength coverage is in- 
sufficient to obtain reliable values for T h ' oZ and L' bol , a flag- 
in the table indicates this and appropriate upper or lower 
limits are given (see the Appendix for more details). 

3.4. Overall Statistics 

Combining all the clouds, we have the following statis- 
tics based on the final c2d data release. The full catalogs 
contain a total of 4.26 x 10 6 entries, of which 6.14 x 10 5 
are also in the high reliability catalog, which requires de- 
tection in at least one band with S/N > 7 and a second 
band, if available at that position, with S/N > 5. A to- 
tal of 6.77 x 10 5 have detections in at least three bands 
and can be further classified. This group provides the 
parent sample from which our YSOc are drawn. The 
great majority are classified as stars (3.32 x 10 5 ) or other 
(2.55 x 10 5 ). The "other" category contains sources that 
do not fit any particular template; the great majority are 
probably galaxies. Of those that remain, 2,965 are galaxy 
candidates and 1,086 are YSO candidates. The statistics 
for individual clouds are given in Table [2j 

The processes described in <j3 . II and H3.2I resulted in the 
list of YSOs used in this paper. There are undoubtedly 
still some contaminants in this list, but we believe that 
they are a small fraction. Overall the ratio of YSOs to 
YSOc is 0.94. The ratio ranges from 0.85 in the Lupus 
clouds to 0.99 in Perseus. In both Perseus and Ophiuchus, 
the fraction is high because of the addition of deeply em- 
bedded and saturated sources not in the original YSOc 
list. 

In the end, there are a total of 1024 YSOs in our sample. 
This number represents an order of magnitude increase 
over previous samples, and they have been selected in a 
uniform way from data with very similar sensitivity. 

4. STAR FORMATION EFFICIENCIES AND RATES 

We discuss these topics first because they rely only on 
the counts of YSOs and the masses of the clouds obtained 
from the extinction maps. As such, the conclusions do not 
depend on further classification or the choice of whether 
or not corrections for extinction are made. 

In Table [H we list the number of YSOs in each cloud, 
along with the number per solid angle and the number 
per area (pc -2 ). The number of YSOs per area is clearly 



highest in Serpens at 14 pc -2 with Ophiuchus second and 
Perseus, Lupus, and Cha II somewhat similar. The lat- 
ter point is a bit surprising because Perseus has many 
more YSOs, but it also covers the largest area (see Ta- 
ble [T] and Figure [TJ). While Cha II is low on YSOs, it also 
has the smallest area. The distribution of YSOs in Perseus 
(|Lai et al.ll2008l ) reveals a large, central section of the cloud 
with almost no YSOs (see Fig. [2]). Together with other 
hints that Perseus may be two overlapping clouds [see dis- 
cussion in lEnoch et al.l (|2006f) ]. this distribution suggests 
that smaller areas for each piece of Perseus could provide 
a more relevant comparison. 

We also estimate the star formation rate from the num- 
ber of YSOs by assuming a mean mass of 0.5 Mq and a pe- 
riod of 2 Myr for star formation. As discussed in |j5l the as- 
sumption of 2 Myr is the estimate of the time taken to pass 
through the Class II SED class, the latest class to which 
our study is reasonably complete. We assume a mean mass 
of 0.5 My, co nsistent with studies of the I MF (IChabrierl 
120031 iKroupal 12001 [Mnkovic fc Traikovskall2006D . There 
may be variations from cloud to cloud , though sma ll num- 
ber statistics are an issue. In Cha II, ISpezzi etafl (pOOl 
derive a mean mass of 0.52 ± 0.11 M Q based on spectro- 
scopic data. The mean s tellar mass may be closer to 0.2 
M Q in the Lupus clouds (jMerin et al.l l2008eJ) . The mean 
mass of the YSOs in Serpens is 0.69 to 0.7 3 M<^, depend- 
ing o n which evolutionary tracks are used (|Qliveira et alj 
[20081) . 

Perseus has the highest star formation rate of 96 
Mq Myr -1 , slightly higher than the rates for Ophiuchus 
and Serpens, while Cha II is conspicuously low at 6.5 
Mq Myr -1 . If normalized by cloud area, however, the dif- 
ferences are less striking (see Table [3]), ranging from 0.65 
M Q Myr -1 pc -2 in Cha II to 3.4 M Q Myr -1 pc -2 in Ser- 
pens. The rate per area in Serpens may be higher in part 
because we covered completely only the part of the cloud 
with Ay > 6 mag, rather than 2 mag for the other clouds. 

The number counts of YSOs do not include any correc- 
tions for unresolved binaries. If a fraction / of all YSOs 
are unresolved binaries, the actual number of forming stars 
and the star formation rates should be multiplied by the 
factor (1 + /), if we ignore even higher multiplic ity. Esti- 
mate s for / range from 0.3 (Lada 2006a) to > 0.5 (M athieul 
11994 . Given that the total stellar mass is always much less 
than the cloud mass, the star formation efficiency would 
be increased by about the same factor. 

4.1. Comparison to Predictions from Kennicutt Relations 

Given the mass surface density of the cloud from our ex- 
tinction maps, one can predict the star formation rate sur- 
face density from the re lations employed for other galaxies 
using the formula from iKennicuta (|1998f l . 

E(SF)(M yr -1 kpc -2 ) = 

(2.5 ± 0.7) x lO -4 (S( 5 as)/lM pc -2 ) 1 - 4±0 - 15 (2) 
Taking all clouds together, the surface density is 64 
Mqpc -2 , and the predicted star formation rate would be 
0.08 M Myr -1 pc -2 , a factor of 20 below the observed 
value of 1.6. As shown in Figure [4] all clouds lie well above 
the prediction of Equation 2, even though these clouds are 
forming only low mass stars, which would be largely invis- 
ible to the kinds of tracers, such as Ha emission, used to 
establish Equation 2. 
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This difference is not surprising, but it reminds us that 
the Kennicutt relation applies to averages over much larger 
regions than individual clouds. The coefficient and expo- 
nent in such relations are clearly functions of the mean 
density of the region being averaged over. A census of 
all star formation within 500 pc should come from the 
completion of many Spitzer observations; a first crude es- 
timate, along with information on the surface density of 
gas, gave rough agreement with t he predictions of Kenni- 
cutt's relation (|Evans et al.ll2008fh N ote that 85% of the 
;as in this region is atomic hydrogen. iBlitz fc Rosolowskvl 
20061) give a more detailed discussion of these points. 

Studies of star formation in dense gas have shown a lin- 
ear relation between star format ion rate and the amount 
of gas traced by HCN emission (iGao fe Solomon I l2004ab 
IGao fc SoWonT 2004b: Wu et al.ll2005ir " As discussed be- 
low, we find that star formation is sharply confined to 
the dense gas, so it is two steps removed from the scales 
on which Equation 2 applies: the formation of molecu- 
lar clouds from atomic gas; and the formatio n of dense 
cores from molecular clouds (|Evans et al.ll2008[ ) . At which 
of these steps does the non-linearity (an exponent that is 
greater than unity) of Equation 2 enter? 

To address that question, studies that resolve individ- 
ual clouds in other galaxies are needed. Studies of M51 
with 0.5 to 2 kpc resolution continue to show a non- 
linea r relation with exponent around 1.4 (|Kennicutt et alj 
2007). Recent observations with sub-kpc resolution of 
both molecular an d atomic gas in other nearby galaxies 
( Big iel et al.1 120081 ) reveal a linear relation between star 
formation rate and molecular gas density over a range 
from 3 to 50 M pc~ 2 . 

£(SF)(M yr^kpc- 2 ) = 

lo -2.1±O.2 (s(jff2)/loM0pc - 2) l.O ± O.2 (3) 

(Note that this relation is normalized to a different 
surface density.) This relation is also shown in Figure 
U] over the range where it was established and extrap- 
olated with a dotted line. Equation 3 would predict 
0.05 1 Mp) yr~ 1 kpc~ 2 , m uch less than our observed value of 
1.6. iBigiel et alj (|2008D note that they are still measuring 
the filling factor of clouds rather than resolving structure 
within molecular clouds. This comparison suggests that 
clouds with the properties we study are filling only about 
3% of their beams. As resolution improves further, and 
as more studies become available within our own galaxy, 
especially of regions forming high mass stars, further com- 
parisons will be important. 

We also plot the relation found for dense g traced 
both in other galaxi es and in massive dense cores in our 
Galaxy by HCN, bv IWu et all (|2005f ). This relation ex- 
trapolated to lower surface densities comes closest to our 
observed po ints. 

Following lElmegreenl (|2002f ). a linear relation in molec- 
ular gas would correspond to a threshold core density for 
star formation of 10 5 cm~ 3 , in good agreement with typi- 
cal densities in the dense cores. However, he also predicts 
that the fraction of total gas in dense cores would be 10~ 4 , 
whereas the fraction in the three clouds with the relevant 
data is 4.6%. Allowing for the fact that 85% of the gas 
in the local kpc is not molecular gets the fraction down 
to 6 x 10~ 3 , still on the high side. The non- linear relation 



that seems to ap ply on larger scales was interpreted by 
lElmegreenl |2002) to mean that the critical criterion is a 
ratio of core density to mean ISM density of 10 5 . Caution 
is still warranted as the extragalactic relations need fur- 
ther analysis, and other factors, such as shear and bars, 
will affect star formation, esp ecially in circum-nuclear re- 
gions (e.g- Uogee et alj|2005| ). 

4.2. Star Formation Efficiency 

Comparison of the mass in YSOs to the cloud mass gives 
a measure of the current efficiency. Since we are only sensi- 
tive to YSOs with infrared excess, this efficiency represents 
an average over the last 2 Myr. We give the star formation 
efficiency, defined as 



SFE = 



M* + M {cloud) ' 



(4) 



in Tabled! The values range from 3% to 6% over the vari- 
ous clouds, with a value of 4.8% taking all clouds together. 
Given the star formation rate (M*) from Tabled we can 
compute a depletion time for the cloud: 

tdep — M (cloud)/ A'h. (5) 

As shown in Table [4] these td ep are 30 to 66 Myr. 
Estimates of cloud lifeti mes range from 10 to 40 Myr 
(Mc Kee fc Os trikcr 2007). If clouds produce stars at the 
current rates for 10 Myr, the final efficiency when star for- 
mation has ended in these clouds could be as high as 15% 
to 30%; as discussed in H4.41 achieving such high efficiences 
would require continued conversion of cloud material into 
dense cores. 

4.3. The Speed of Star Formation 

iKrumholz fc T an ( 2007| have emphasized that star for- 
mation, even in molecular clouds, is slow in the sense 
that the star format i on rat e per unit free fall time is low. 
Krumholz & McKec (2005) define the star formation rate 
per unit free fall time (SFRff) to be the fraction of an 
object's mass that turns into stars in a free fall time at the 
object's mean density. Using our cloud masses from the ex- 
tinction maps covering the same areas where we surveyed 
for YSOs, we can quantify this, again where we mean star 
formation over the last 2 Myr since we identify only the 
sources with infrared excess. Then we have 



SFR ff = M+t f f/M(cloud) = t ff /t de 



(6) 



where tff is the free fall time for the mean density of the 
cloud, calculated from 

tff = 34Myr/7(^). (7) 

The density (n) in that calculation is the density of all 
particles, and we have assumed a mean molecular mass of 
2.3 amu. Mean densities for the clouds, computed from the 
mass and surface area by assuming spherical clouds, are 
given in Table [TJ The average over all clouds is (n) = 390 
cm' 3 . Then, we find SFRff ranging from 0.028 to 0.064, 
with an average over all clouds of 0.040. This value is on 
the high side of the ran ge of SFRff inferred fr om more 
global co nsiderations by IKrumholz fc Tanl (|2007l ) and the 
theory of IKrumholz fc McKed (|2005h . 
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4.4. Efficiencies and Rates in Dense Gas 

Figure shows a strong correlation of YSOs, especially 
the younger ones, with regions of high extinction. For the 
clouds with millimeter continuum maps, a strong corre- 
lation of the youngest objects with d ense cores identifie d 
in the continuu m maps is apparen t (jEnoch et al.l l2007t) . 
As discussed bv lEnoch et al.1 (|2007j ). the cores identifed in 
this way have minimum mean densities of about 2 x 10 4 
cm~ 3 and more typically more than 10 5 c m" 3 , and al- 



most a ll seem to be gravitationally bound ([Enoch et alj 
I2008bh . Their mean densities are 50 to 200 times the 
mean density of the cloud; they thus represent quite dis- 
tinct entities, rather than just modest peaks in column 
density. Thus, it is also interesting to calculate the ef- 
ficiencies and rates using the total mass of gas in dense 
cores for the three clouds wi th Bolocam data, using the 
data from lEnoch et all (|2007t) . The total mass in YSOs in 
the cloud is quite comparable to the mass in dense cores 
(Table [4]); taking all three clouds together the ratio is 1.3. 
We compute the time to deplete the current stock of mass 
in discrete dense cores by dividing the SFR by the mass in 
dense cores. We find 0.6 Myr for Ophiuchus, 1.6 Myr for 
Serpens, and 2.9 Myr for Perseus. Taking all clouds to- 
gether, tdep(dense) — 1.8 Myr. These times are consistent 
with the 2 Myr timespan for detectability of YSOs that 
we use to calculate lifetimes and with plausible spreads of 
formation times in clusters. While star formation is faster 
and more efficient in the dense gas probed by Bolocam, it 
is still slow compared to a free fall time. Using a mean 
density of 5 x 10 4 cm -3 to calculate the free fall time for 
the typical core leads to SFRff (dense) — 0.05 for Perseus 
to 0.25 for Ophiuchus. 

These values for SFRff (dense) are somewhat higher 
than SFRff for the whole cloud, but they neglect the 
fact that the mass of stars is now comparable to the mass 
of dense gas. The depletion times also assume that all the 
gas winds up in st ars once it r e aches the density of the 
cores identified by lEnoch et al.l (|2007f ) . We can make a 
second calculation in which we assume that star forma- 
tion began at some time in the past with a larger reservoir 
of dense gas and that a fraction e of the dense core mass 
winds up in a star, while a fraction 1 — e is lost to the star 
forming region by outflows. By comparing th e core mass 
functi on to the initial mass function of stars, Alves et alj 
(|2007[) estimated that e = 0.3 and lEnoch et aU (|2008bD 
finds e > 0.25. At any given time, tdep — M (dense) /M* 
and the mass of dense gas decreases exponentially, with a 
time constant of etd ep - At time t after the start of star 
formation, the expression for tdep becomes, 



eln(l + rj/e) 



(8) 



where n = M^/M (dense) at time t. As usual, we take 
t = 2 Myr as the timescale over which we have complete 
statistics, assume e = 0.3, and use the observed value of 
n. This expression then yields values for tdep that are 
more like 3-6 Myr and SFRf f (dense) that are 0.03 to 
0.06, comparable to th ose for the cloud as a whole, as 
iKrumholz fe Tanl (|2007t ) would predict. 

5. CLASSES AND LIFETIMES: THE STANDARD ANALYSIS 



5.1. Historical Classification Methods 

The current working mo del for star formati on arose in 
the 1980s with the merger (Ad ams eTaII[l987h of an em- 
pirical classification scheme, based on the slope of the 
spectral energy distrib ution (SED) between 2 and 20 /an 
(jLada fc Wilkindll984h . with a theoretical picture of star 
formation i nvolving the collaps e of an isolated rotating 
dense core (iTerebey et al.lll984f) forming a star and disk 
(Ad ams fc Shul l 1986). The essential stages of the theoreti- 
cal model were three: the collapse of the envelope, forming 
the protostar and disk; the continued accretion of disk ma- 
terial onto the forming star; and the dissipation of the disk 
by planet formation, evaporation, etc. Each of these stages 
became identified with an empirical SED class as follows: 
the envelope collapse as Class I; the accretion disk and star 
as Class II; and the dissipation of the disk during Class 
III. With further study, other significant events could be 
distinguished theoretically, and further refinements of the 
classifi cation system we r e sugg ested. We follow the sugges- 
tion of iRobitaille et all (|2006f ) in referring to the physical 
arrangement as a "Stage" and the SED characteristic as a 
"Class". 

We focus now on the evolution of the e mpirical classifica- 
tion sy stem. Following the recognition bv lLada"fc Wilking 



1984) that SEDs were falling into natural groups, lLada 



19871) first codified the tripartite class system using the 



spectral index. The original boundaries were as follows: 

I < a<3 

II -2<a < 

III -3 < a< -2 

Note that Lada and coworkers consistently use the sym- 
bol a, but we have standardized on a. These were sup- 
plemented by descriptions; for example, Class III sources 
could have some "mid-infrared excess" , but "no or little ex- 
cess near-infrared emission." Class II and III objects were 
all visible, but Class I objects were all invisible for A < 1 
/im. Class II sources were identified as T Tauri stars. It 
is worth noting that this system and many of the later 
elaborations were based on study of the cluster of sources 
in the L168 8 (Ophiuchus) c l oud. 

By 1989, IWilking et ail (|1989h had noted that some 
SEDs had "double humps", one around 1 fj,m and one 
in the far-infrared. These were put into subclasses, la- 
beled II-D and III-D. Class II-D sources were T Tauri 
stars with far-infrared excesses and were thought to be 
intermediate between Classes I and II. Class III-D sources 
looke d like redde ned photospheres plus a far-infrared ex- 
cess (lLadalll991l) . By comparing the n umbers of objects 
in different classes. IWilking et~a l. (1989) also assigned rel- 
ative lifetimes to the different classes (discussed below). 

The next major develo pment was the intro duction of a 
class "before" Class I by I Andre etafl (fl99l . Since zero 
was not part of the original Roman number system, the 
Arabic symbol for zero was used. Since the Class objects 
could not then be obs erved at the wavel engths originally 
used for classification, lAndre et al.l (|2000D listed three cri- 
teria for a Class source: 

1. Indirect evidence for a central YSO, as indicated 
by, e.g., the detection of a compact centimeter ra- 
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dio continuum source, a collimated CO outflow, or 
an internal heating source; 

2. Centrally peaked but extended submillimeter 
continuum emission tracing the presence of a 
spheroidal circumstellar dust envelope (as opposed 
to just a disk); 

3. High ratio of submillimeter to bolometric luminos- 
ity, suggesting that the envelope mass exceeds the 
central stellar mass: L smm / 'Lboi > 0.5%, where 
Lsmm is measured longward of 350 /im. In prac- 
tice, this often means an SED resembling a single- 
temperature blackbody at T ~ 15 — 30K. 

The boundary between Class and Class I was asso- 
ciated with a new distinction between physical stages 
(jAndre et al.lll993T) : the point at which the masses of the 
protostar and the remaining envelope were about equal. 

The last m ajor development came in 1994, when 
IGreene et al.l (|1994f ) formalized the 4-class svstemFI. 
again based on the L1688 cluster, as follows: 

I 0.3 < a 

Flat -0.3 < a < 0.3 

II -1.6 < a < -0.3 

III a < -1.6 



mean frequency as the actual SE D, considering all wave- 
lengths with available data. I Chen et al.l (|l995f ) showed 
that the traditional classes, including 0, could be associ- 
ated with certain ranges of T&oZ, as follows. 

T bo i < 70 



IGreene et al.l (|1994| ) described the evolutionary status 
of sources in the new Flat Class as "uncertain", but 
ICalvet et al.l (| 19941 ) s howed that these cou ld be interpreted 
as infalling envelopes. IGreene et al.l (|l994h noted that they 
found no sources with —1.9 < a < —1.6 or —0.3 < a < 0. 
These classifications were based on the slope between 2 and 
10 fini, though they claimed that comparison to classes 
based on slope s between 2 and 20 /zm showed that this 
did not matter. IGreene et al.l d~1994) did not list Class in 



their list, but they did acknowledge lAndre fc Montmerld 
(|1994h . who pointed out that millimeter wavelength emis- 
sion became much weaker for a < —1.5, as a rationale for 
the revised boundary between Classes II and III. 

Improving sensitivity at millimeter wavelengths 
led to the detection of some starless cores 
(| Ward-Thompson et al.|[l994T ). and these were called "pre- 
protostellar cores" (PPCs) or later "prestellar cores". To 
fit t hem into the class s ystem, one might use Class — 1, 
but iBoss fc Yorkd (|l995l ) argued that they should instead 
be Class —2 to make room for a theoretically important 
event, the formation of the first (molecular) hydrostatic 
core, which would then be called Class —1. These terms 
have not caught o n , and the recent definitive review by 
Idi Francesco et all (|2007l ) simply distinguishes "prestel- 
lar" cores among the larger set of starless cores as being 
gravitationally bound. 

With the increasingly baroque nature of the class sys- 
tem, it was natural to seek a continuous variable to capture 
the transitions. By analogy wi th the effective t emper ature 
of stars and spectral classes, iMvers fc Laddl (|l993l ) sug- 
gested the use of a bolometric temperature, defined as the 
temperature of a blackbody with the same flux-weighted 

17 The class of sources with a = 0.3 was left undefined in Greene et 



I 70 < T bo i < 650 

II 650 < T boi < 2800 

For very early stages, data at short wavelengths were 
lacking, and Tboi could be quite sensitive to whether or 
not such data existed. Detections at shorter wavelengths 
could increase Tboi enough to move a source from Class I to 
II. Orientation of aspherical envelopes and disks can also 
affect Tboi substantially. At the other wavelength extreme, 
improved availability of submillimeter data considerably 
increased the number of Class sources (|Visser et al.|[2 002. 
lYoung eraT1l2003h . 

For distinguishing Class from Class I, the ratio of 
luminosities [Lboil Lsmm) may be most useful. In this 
scheme, a ra t io exc eeding 200, the inverse of criterion 3 of 
I Andre et all (2000). marks the transition from to I. This 
form of the ratio has the virtue of increasi ng with the evo- 
lution ary progression, as does Tboi , and lYoung fc EvansI 
(2005) found that this ratio was a more robust indica- 
tor than Tboi of the ratio of mass in the star to mass in 
the envelope , based on models of colla psing, cores. Both 
IVisser et all ([20021) and lYoung et all (|2003j) found that 
quite a few objects would be classified as by the ra- 
tio, but as I by Tboi- A drawback to using Lboil L smm is 
the lack of complete data at 350 /im, the shortest wave- 
length included in L smm ; we lack such data for many of 
our the sources in this study. Clearly, conclusions about 
evolution will depend on establishing a clear connection of 
such parameters to physically relevant stages of evolution. 

It remains difficult to capture the increasing amount of 
information in any single parameter. With the greatly 
increased and more uniform data set available from the 
c2d project, we will reevaluate the different tracers to 
learn what, if any, well defined criteria can trace evolu- 
tion through all stages. The distinction between physi- 
cal stages of evolution and SED classes, however defined, 
has be en us efully emphasize d recen tly by iRobitaille et alj 
(|2006h and IRobitaille et all (|2007l ). who presented large 
grids of models SEDs from 2-D r adiative transfe r calcu - 
lations, following earlier work bv IWhitnev et all ([20031) . 
which de monstrated the impo rtance of inclination effects. 
Similarly, ICrapsi et al.l (|2008l ) have discussed the possibil- 
ity of confusion between stages and classes in the context 
of a grid of 2-D models. 

5.2. Lifetimes: Previous Estimates 

One of the goals of research into star formation is to 
constrain the lifetimes ass ociated with di f ferent stages of 
evolution. As noted above. IWilking et alj (|1989| ) used the 
numbers of objects in various classes to establish relative 
lifetimes. This method assumes that the census is com- 
plete and that star formation in the sample has been a 
continuous process at a steady rate for longer than the 

ID {Till), so we have arbitrarily assigned such sources to Class I. 
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dwell time in the slowest phase considered. When extrap- 
olated from a single region, it also assumes that other vari- 
ables are irrelevant, such as the total mass available, the 
conditions in the region, such as turbulence, etc. 

With these caveats in mind, IWilking et al.1 (|1989t ) used 
their data on the Ophiuchus cluster to suggest lifetimes. 
They estimated the ages of the Class II sources, finding an 
average of 0.39 ± 0.17 Myr, with the oldest at about 1.5 
Myr. They found roughly equal numbers of Class I and 
Class II sources, suggesting equal lifetimes. With some 
other considerations, they suggested a duration for the 
Class I phase of 0.2 to 0.4 Myr. Later. iGreene et all (|l994h 
found that Class I plus Flat SED lifetimes in Ophiuchus 
were 75% of the Class II lifetime, which they took to be 
0.4 Myr. A somewhat older population in O phiuchus, with 
ages a round 2 Myr, was later established bv lWilking et alj 
(2005^ 



Kenvon fc Hartmannl (|l995l ) analyzed data from the 



Taurus cloud and found about ten times as many Class 
II plus Class III sources as Class I sources, very different 
from the situation in Ophiuchus. They found a similar 
number of Flat SED sources and suggested lifetimes of 0.1 
to 0.2 Myr for each of the Class I and Flat classes, based 
on a duration for Class II plus Class III of 1 to 2 Myr, 
which they found from comparison to evolutionary tracks. 

More recent analyses of the Perseus cloud find a lifetime 
for the Clas s 0/1 SEP of 0.25 to .67 Myr within 95% confi- 
dence level (jHatchell et al.|[2007t ). These estimates include 
Spitzer data on embedded sources a nd the analysi s of th e 
age spread of the IC348 cluster from lMuench et ail (|2007f ). 
Early estimates of 0.01 Myr for the lifetime of the Class 
phase were bas ed on the small nu mber of Class sources 
in Ophiuchus (|Andre et alj Il99"3l) . Because the first ex- 
amples of Class sources found were quite luminous with 
powerful, collimated outflows, the short lifetime suggested 
a phase of rapid accretion i n which at least ha lf of the final 
stellar mass was accreted ([Andre et al.ll2000l ). 

Most of these estimates were based on samples of 50 to 
100 sources and in particular clouds. They are subject to 
small number statistics and to possible differences from 
cloud to cloud, especially for the earlier classes. An excep- 
tion was a compilation of 95 Class 0/1 objects (|Froebrichl 
2005). The positions of these sourc es in Th n i-Lh„i plane was 
compared to grids of models, and iFroebrich et al.1 ((2006) 
favored a Class lifetime of 0.02 to 0.06 Myr. 

5.3. Numbers and Lifetimes: The c2d Large Cloud 
Sample 

Table [5] present s the numbers of YSOs in each (2008) 
I Greene et all (|l994h class as determined by 2MASS - BothQ 
Spitzer 24 /zm photometry. They are also s hown graphi- 
cally i n Figure [5l We remind the reader that lGreene et all 
(|1994h used photometry only out to 10 /urn, but found no 
significant differences in a for sources with data out to 20 
/xm. We have not yet separated Class sources from Class 
I sources in Table [5] this separation will be discussed in 
H7.H In addition, Class III is certainly missing sources, as 
discussed in Sj3j So at this point, we restrict comparison 
to Class I, Flat, and Class II sources. 

With 1024 YSOs and more than 100 in each class, the 
uncertainties from small number statistics are decreased to 
less than a 10% effect in any class, almost certainly smaller 



than other sources of uncertainty. We find that 60% of our 
YSOs are in Class II, with 16% in Class I, 12% in Flat, 
and 12% in Class III. The results with classes defined by 
a' (after extinction corrections) are also shown in Table 
[5] for each cloud and summarized in Table [H The main 
effect of extinction corrections is to decrease the number 
of Class I (to 14%) and Flat (to 9%), while increasing the 
number of Class II (to 64%) and Class III sources (but still 
at 13%). 

With or without e xtinction corre ct ions, our result is 
quite unlike that of IWilking et al.1 (fl989) who found 
roughly equal numbers of Class I and Class II objects in 
Ophiuch us. We find m a ny m ore Class II objects, as also 
found by Wi lking et al.l (|2005h . The variations from cloud 
to cloud are substantial, as can be judged from Figure 
though the effects of small number statistics are substan- 
tial for Cha II and Lupus (see Table [5]). Among the three 
clouds with substantial numbers of sources in early classes, 
Perseus stands out as particularly rich in Class I sources. 

If star formation has been continuous over a period 
longer than the age of Class II sources and if we can av- 
erage all our clouds, we can obtain relative lifetime es- 
timates for each phase by taking the ratios of number 
counts in each class, and multiplying by the lifetime for 
Class II. Recall that this method assumes that all objects 
of all masses behave identically and that no other variables 
enter importantly. Since these assumptions are not very 
realistic, caution is needed in interpreting these, or any 
previous, estimates for lifetimes. 

The lifetime of the Class II phase i s still uncertain. In 
the classic study of 83 cTTS in Taurus. IStrom et all (|l989h 
found that 60% of those with ages less than 3 Myr had 
K-band excesses. The sample of cTTS would probably 
be bia sed towards those w ith long-lived disks. More re- 
cently, Hais ch et al.l (|200lh studied young clusters of dif- 
ferent ages and showed that half the stars had lost their 
disks, as indicated by an L-ban d excess, in < 3 M yr. A 
study of disk frequency in wTTs (|Cieza et alj|2007t ) shows 
that half have lost their disks within about 1 Myr; since 
the wTTS sample is undoubtedly biased towards stars 
that lose disks early, this is probably a lower limit for 
the disk half-life. The median age of stars in IC348, as- 
suming constant form ation rate and the d istance we have 
adopted, is 3.0 Myr (|Muench et al.ll2007f ). The ratio of 
wTT S to cTTS in IC34 8 is 1.5, a nd somewhere between 
30% (|Lada et al.ll2"006bh and 22% (|Cieza fc Baliberl[200l 
of the wTTs have disks. However many of these would not 
be Class II objects. These numbers suggest a n upper limit 
of 3 M yr for the Class II phase. For Cha II, ISpezzi et al.l 
find a mean age for the YSOs of about 2 Myr. 
Both iBontemps et al.1 (|200lD and iKaas et alj (|2004h also 
assumed 2 Myr in their analysis of ISO data. Based on all 
these pieces of information, it seems that the best estimate 
for the lifetime of the Class II phase is 2 ± 1 Myr. The 
uncertainty is dominated by uncertainties in stellar ages, 
whic h are uncerta i n by fa ctors of 2 at these early times 
[e.g., lHaisch et all (|2001l ). iHillenbrandl (|2008l )]. Since the 
marker used to define the transition was that half the stars 
lacked infrared excess, this lifetime may be best thought 
of as a half-life rather t han a lifetim e which is the same 
for all objects [see also IHillenbrandl 



In partic- 
ular, there is some evidence for longer lifetimes for in- 
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frared excesses aro und lower mass stars or brown dwarfs 
(| Alters et alJl2007l) . 

If we take 2 Myr for the duration of the Class II phase, 
then the lifetime of the Class I phase, t(I) =0.54 Myr, with 
an additional t(F) =0.40 Myr for the Flat SED phase. 
Using the numbers after extinction corrections, we obtain 
t'(I) =0.44 Myr and t'(F) =0.35 Myr. With or with- 
out extinction corrections, these estimates are substan- 
tially longer than som e recent estimates f or the Class I 
lifetime of 0.1-0.2 Myr (|Greene et al.lll994h . but this esti- 
mate depends on the assumption of an age for the Class 
II phase. Interestingly, our estimate agrees better with 
the original es t imate of 0.39 Myr for the Class I stage by 
IWilking et al.l (|l989f ) because they assumed an age of 0.4 
Myr for the Class II stage in the L1688 cluster and found 
equal numbers. So we obtain the same answer but with 
totally different data and assumptions. 

The uncertainties in these numbers are clearly not dom- 
inated by small number statistics any longer. While the 
original definition of classes did not include extinction cor- 
rections, it appears that an approximate correction leads 
to 10% and 20% shorter lifetimes for the Flat and Class 
I phases. The limitations of the underlying assumptions 
dominate the uncertainties. In particular, the continuous 
formation assumption is dubious. If star formation in a 
particular region decreased or stopped within the last 2 
Myr, we would underestimate the lifetime for the Class I 
and Flat classes. Conversely, if star formation began or in- 
creased within the last 2 Myr, we would overestimate those 
lifetimes. Evidence for these effects can be seen in the life- 
times calculated cloud by cloud in Table 03 For Class I, 
these vary from t(I) = 0.2 Myr (Cha II and Lupus) to 
t(I) = 0.77 Myr in Perseus, which supplies over half the 
Class I objects and nearly 40% of all YSOs in our sam- 
ple. Since Perseus clearly had significant star formation in 
the IC348 region 3 Myr ago, it does not obviously violate 
the continuous formation assumption, but there may have 
been a recent burst. If we take an extreme position and 
leave Perseus out of the average, (t(I)) = 0.40 Myr. 

Our only defense against these issues is to hope that 
averaging over all the clouds provides some cancellation 
of biases. Since we selected our clouds to have a range of 
known star formation activity, the sample should not be 
too heavily biased, but it does contain three of the most ac- 
tive nearby regions. Once the data from the Taurus cloud 
and the Gould Belt Survey are available in the same form 
as the c2d data, these numbers should be recomputed. 

One possible systematic bias towards overestimating 
Class I lifetimes would arise if low mass objects are found 
in our sample when they are young and accreting but 
fall into the galaxy confusion region during the Class II 
stage. We start to lose Class II sources around the stel- 
lar/substellar boundary, based on a few examples. If we 
assume that we miss all substellar objects, which is too ex- 
treme, we would miss 10% to 30% of Class II sources (e.g., 
lAndersen et~ail 120061 . iLuhman et aLll2007t ). which would 
introduce a comparable percentage error in the Class I 
lifetime. 

6. ALTERNATIVE CLASSIFICATIONS 

In this section, we raise s ome questions abou t the use of 
the class system defined bv lGreene et al.l (|1994|) to catego- 



rize YSOs. First we retain a as the discriminant, and dis- 
cuss uncertainties attendant on the choice of wavelengths 
and method of calculating a. We then discuss alternatives 
for very early and very late stages. 

6.1. Definition of Classes with a 

The original definition (|Lada fc Wilkind I1984D used 
wavelengths between 2 and 20 fim to define a. We use 
a least-squared fit to a ny data in Table IA8I between 2 and 
24 /Ltm. More recently, lLada et "all (|2006bl ) have used a fit 
to only the 4 IRAC bands (3.6 to 8.0 /im) in a study of 
disks in IC348. How much difference does the choice of 
wavelength range make? We classified sources both ways 
(Tabled . Using only IRAC bands moves sources from ear- 
lier to later classes in such a way that we get more Class 
III (by 47%), fewer Class I (by 8%), fewer Flat SEDs (by 
15%), and fewer Class II (by 8%). Reddening will have 
a larger effect at shorter wavelengths, which biases a up- 
ward, producing "Flat" SEDs that are clearly reddened 
Class II sources when the full SED is examined. Thus, dif- 
ferent definitions of a can result in 10% to 15% difference 
in inferred lifetimes for the earlier SED phases. We did 
another experiment, using sources in C ha II, for which we 
have spectral types (jSpezzi et al.ll2008l ). allowing photom- 
etry to be corrected for extinction. Indeed, the values of a 
calculated after extinction corrections were smaller, with 
a mean difference of 0.2. This was not enough to affect 
the classification of sources in Cha II. With the approxi- 
mate extinction corrections used for the full data set, we 
find somewhat larger effects, but still 10% to 20%. These 
effects are probably smaller than other sources of uncer- 
tainty, inherent in assumptions such as the continuous for- 
mation assumption f £|5.2[) . 

6.2. Definition of Classes with T^oi 

Values for L^oi and T^oi were calculated as described in 
H3.3I and are given in Table IA12I The observed values 
of Lboi and Tf, i are easily compared to most other data 
on embedded sources in the literature, but i n the paper 
conne cting T^oi to the classes defined by a, IChen et al.l 
( 1995]) corrected the observed flux densities for extinction 
before computing T(, i- For less embedded objects with 
Ay determinations from optical data, they used those val- 
ues. If a source had no measured Ay but had a known 
spectral type, they used the spectral type and V — R to 
calculate Ay. If no spectral type was known, they as- 
sumed an M0 spectral type. For sources with no optical 
data, they assumed an average Ay — 1 for Taurus and off 
the core region in Ophiuchus, and Ay = 10 for the core 
sample in Ophiuchus. Our extinction corrections are sim- 
ilar to what they assumed in Ophiuchus and substantially 
larger in other clouds (see £13. 3p . They analyzed the ef- 
fects of extinction; these could be substantial in both T^oi 
and Lb i, primarily for the sources with higher intrinsic 
Tboi, where most of the energy emerged at shorter wave- 
lengths, and these sources generally had reasonably well 
determined values for Ay. They considered the uncer- 
tainties introduced by extinction uncertainties and found 
a factor of two for high T^ oi sources and 10% for low Tb i 
sources (IChen et alJll995li . Most subsequent application 
of Tboi has been to deeply embedded sources, so no cor- 
rections for extinction are typically made; the energy ab- 
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sorbed at short wavelengths is assumed to be re-radiated 
at longer wavelengths. 

Our method of exti nction corrections attempts to cap- 
ture the essence of the iChen et all (fl995h approach. If we 
compare the numbers in classes defined by T' hol , the effects 
are more dramatic than was the case for classes defined by 
a. The observed fluxes produce only one Class III source 
with the Tboi dividing line from lChen et alj (|l995f ). while 
using T' hol give 485 Class III sources! Clearly, the bound- 
ary in T{, i between Class II and III is a sensitive function 
of extinction. 

Class objects cannot be separated from Class I objects 
by using a. Previously this statement was true because 
Class objects cou ld not be detected in the mid-infrared 
(| Andre et al.ll2000D . With IRAC on Spitzer, we now de- 
tect many, though not all (jJprgensen et al.| [2005). Class 
sources. However, the resulting values of a from the c2d 
catalog scatter widely and do not correlate w ith other dis- 
crimin ants such as Tboi for the early classes ([Enoch et alj 
2008a|). We show the full sample of sources with enough 
data to constrain both a and Tboi in Figure [6] While 
there is some tendency for Class sources, as identified 
by Tboi, to have larger values for a than Class I sources, 
the scatter precludes any Class criterion based only on 
a. These conclusions do not change using the photometry 
after extinction corrections to calculate Tboi (right panel 
of figure). The range of a may be due in part to geometry 
of th e envelope-disk interfa ce, the inner radius of the enve- 
lope ( J0rgensen et al.ll2005h . different amounts of scattered 
light, emission fro m jets in the 4.5 /im band, and deep ice 
features (see, e.g..|Huard et alj|2006t iDunham et al.| [2006: 
iBourke et al.ll2006l iBoogert et al.ll2008D . The fluxes do not 
increase monotonically from 3.6 to 24 ^m in some Class 
sources. 

Conversely, there was no definition of the Flat SED 
class in terms of Tboi ■ Figure [5] shows a strong concen- 
tration around 65 K, the Class I/II boundary defined by 
IChen et alJ (fl99l . with some outliers. Histograms of the 
number of sources in linear bins of Tboi and T 6 ' oi are shown 
in Figure [7] The different classes according to a and a' are 
color-coded. The mean value of Tboi for Flat SED sources 
is 649 K. We suggest a range of Tboi from 350 to 950 K 
for Flat SED sources; this range includes 79% of the Flat 
SED sources and contains only 23% of the C lass I sources 
(by a) and 22% of the Class II /III sources. lEnoch et al.l 
(20083) have also suggested a boundary between Stage I 
and Stage II sources in the 400 to 500 K range. 

The distribution of Flat SED sources in extinction- 
corrected T h ' oi is somewhat broader and shifted to higher 
T' hol , with a mean value of 844 K (right panel of Fig. [S]). 
A range of 500 to 1450 K contains 77% of the Flat SED 
sources, while only 14% of Class I and 18% of Class III 
sources lie in this range. Figure [7] suggests that Flat SED 
sources do form a distinct class of objects. However, this 
seeming homogeneity can in part be caused by observing 
a flat SED over a restricted wavelength range; Tboi has 
little freedom to vary without data at shorter and longer 
wavelengths. 

Class II and III sources have barely distinct peaks in Fig- 
ure but substantial overlap. This figure illustrates the 
sensitivity of the Class II/III boundary to extinction cor- 
rections that was noted above. It also suggests that Tboi is 



not particularly good at distinguishing Class II from Class 
III objects and a is preferred. 

6.3. Source Types from c2d Catalogs 

We can also compare the source types in the c2d cat- 
alogs, as described in fJ51 to classifications based on a or 
Tboi ■ Figure \E\ shows the source types versus Tboi , both 
with and without extinction corrections, with the color 
code showing the a SED class. The categories such as 
YSOc_red, red, and rising are almost exclusively Class 0/1 
or Flat. The source type YSOcj3tar+dust(IRl) contains 
a broad range of classes; an excess starting at 3.6 /xm 
can continue to increase (Class I), stay constant (Flat), 
or fall (Class II or Class III). When the excess begins 
at longer wavelengths, such as YSOc_star+dust(MPl), al- 
most all are Class III, with a few exceptions. The excep- 
tio ns may include the "cold disks," such as those studied 
bv lBrown et al.l {2007). The generic YSOc class contains a 
heterogeneous sample of classes. Similar patterns are seen 
in Figure [51 which shows the source type versus a. 

6.4. Source Types by Color-Color Diagrams 

C olor-color diagram s constructed from IRAC photome- 
try (|Allen et al.ll2004[) or from combined photometry from 
IRAC and MIPS (jMuzerolle et al1l2004l ) have also proved 
useful in separating sources in different phases of evolu- 
tion. These are shown for our YSO data in Figures [TU1 
andfTTl both with data after extinction corrections in the 
right panel. The color codes correspond to the classifi- 
cation by a except that Class and I sources are sep- 
arated by Tboi- The overall correspondence is good, with 
Class II sources heavily concentrated in the expected loca- 
tions in_Hie_twodia^rarns, as indicated by the boxes , taken 
from I Allen et al.1 (|2004D and IMuzerolle et al.l (120041) . F lat 
SED sources we r e not separated by lAllen et al.l (|2004h or 
IMuzerolle et al.1 (|2004h . but they ten d to lie between th e 
Class II and Class 0/1 regions [see also I Allen et al.1 ([2007ft ]. 

The Class /I sources scatter more widely but generally 
have redder colors than the Flat SED sources. The Class 
and I sources are not clearly separated in these diagrams, 
similar to what we found in Figure^ but on average, Class 
sources have redder colors. There are exceptions, such 
as Class sources with [5.8] — [8.0] colors near zero. Also, 
there are some Class II sources rather far from their home 
box, especially in [3.6] — [4.5], many of which find their way 
home in the figure corrected for extinction. 

Our large sample suggests that some adjustments to 
the box locations and perhaps shapes would be advisable 
for better segregation of classes in the [3.6] — [5.8] versus 
[8.0] — [24] diagram. In particular, the divisions between 
Class II, Flat, and Class 0/1 may be captured better by 
diagonal lines in Figure QTJ The bottom two panels of Fig- 
ure QT] show the divisions b etween the simulated S EDs for 
sources in different stages ([Robitaille et al.l [2006h . which 
do capture the divisions seen in the data well, especially 
for the dereddened data. This general congruence indi- 
cates that class definitions generally agree with the phys- 
ical stages, despite exceptions. We consider some of the 
remaining issues and extend the analysis to classes not 
separated in Figure [TT1 (prestellar. 0, and Flat) in the next 
section. 
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7. CONNECTION TO PHYSICAL STAGES 

In this section, we discuss the relations between SED 
classes, defined above and physically meaningful stages in 
the evolution. We will estimate lifetimes for Class and 
prestellar phases. 

7.1. The Early Years: Class and Before 

As discussed previously, distinguishing Class from 
Class I sources can be ambiguous without very well sam- 
pled SEDs through the submillimeter region. Lacking data 
at 350 /im for many of our sour ces, we will use Tp i = 70 
K to define the Class boundary (|Chen et al.lll995ft . 

To include still earlier phases, prestellar cores, we rely 
on surveys in the millimeter continuum. We have complete 
surveys at 1.1 mm with reasonably uniform sensitivity for 
Perseus, Serpens, and Ophiuchus and a less sensitive sur- 
vey of C ha II. Focusing on the first 3 clouds, lEnoch et al.l 
(2008b) report on a total of 201 cores (defined as sources 
of millimeter continuum emission above so me level). By 
compar ing these cores to the Spitzer data, lEnoch et alj 
( 2008a) separate these into cores with and without central 
luminosity sources. The latter are assigned to prestellar 
cores. Cores with central luminosity sources are divided 
into Class 0, I, and II using the b oundaries in Ttoi listed in 
H5.ll based on I Chen et"al (|1995ft . There are only a hand- 
ful of Class II sources associated with dense cores, and 
these are likely due to projection effects. The results are 
given in Table [HI 

To get lifetimes for these classes, we must bootstrap, 
using the lifetime derived above for Class I (now reinter- 
preted as the lifetime for Class and Class I combined). 
Using the same continuous formation assumptions and a 
lifetime for Class 0/1 of 0.54 Myr, we would infer a mean 
lifetime for Class of 0.16 Myr for all three clouds com- 
bined. If we use the same sample, but use T' holl calcu- 
lated after extinction corrections, the resulting t'(0) = 
0.10 Myr. The 60% decrease in lifetime results partly 
from the shorter lifetime for Class I, but primarily from 
the smaller number of Class sources after extinction cor- 
rections (Table [6]). Doing the calculation cloud by cloud, 
with the values of t[T) for that cloud, yields a wide range, 
from Ophiuchus with t(Q) = 0.043 Myr to Perseus with 
0. 32 Myr. The Perseus result is consistent with the results 
of lHatchell et alj (|2007t ). While the Ophiuch us lifetime is 
3 time s higher than the original estimate of lAndre et alj 
41993ft . it is much shorter than for the other clouds. In ad- 
dition, it could be overestimated, as it rests on three Class 
sources, two of which are borderline Class 0/1 sources. 
Either Class sources in Ophiuchus do evolve very quickly 
into Class I sources, or the continuous formation approx- 
imation breaks down there. A decrease in star formation 
rate beginning 0.16 Myr ago would result in very few Class 
sources now. The variation from cloud to cloud is less if 
the s ame value of t(I) of 0.54 Myr is used for all clouds, 
as in lEnocfTet al. (2008a!). Note that these results use the 
ratio of Class to Class I sources with both required to 
have millimeter emission. 

To get a lifetime for the pres tellar (t(pre)) cores th at are 
seen in the Bolocam surveys, lEnoch et al.l (|2008bft com- 
pared the number of cores that remain starless to the 
number of embedded protostars. The latter number is 
slightly less than the number of Class 0/1 sources because 



a small number of Class I objects do not have distinct 
Bolocam cores. Since this number is small, we still use 
the lifetime for the combined Class and I phase to rep- 
resent the lifetime of such protostellar cores. The result, 
using all cores in all clouds, is t(pre) = 0.46 Myr. In 
this case, extinction corrections remove enough Class 0/1 
sources to balance the decrease in Class I lifetime, so the 
extinction-corrected result is t'(pre) =0.47 Myr. These 
nu mbers differ very sli ghtly from that (0.45 Myr) given 
by lEnoch et al.l (|2008bft because of small differences in the 
source list. One must immediately caution that this is 
the lifetime for cores remaining starless after they become 
dense enough to be detected with Bolocam. B ased on care- 
ful analysis of the instrumental limitations, Enoch et al. 
( 2007) concluded that the minimum mean density is about 
2 x 10 4 cm~ 3 . The variation in t(pre) from cloud to cloud 
is 0.23 (Serpens) to 0.78 (Perseus) Myr if we use the t(I) 
for each cloud; Ophiuchus is not unusual in this category 
with a value (0.37 Myr) near the mean. However the total 
mass of the cores in Ophiuchus is small; it may be a cloud 
that has just finished a burst of star formation that has 
exhausted much of its dense gas. Again, the range is less 
if the s ame t(I) is used for all clouds, as in lEnoch et al.l 
(|2008bft . 

Comparison to a separate study using SCUBA data can 
provide some id e a of t he uncertainties in these estimates. 
|j0rgensen et al.l (|2007ft co mpared SCUBA data i n Perseus 
to c2d infrared data, and |j0rgensen et al] (|2008l ) have ex- 
tended this study to Ophiuchus. Compared to the Bolo- 
cam data, the SCUBA data are taken at a slightly shorter 
wavelength (0.85 mm versus 1.1 mm), have a smaller beam 
(15" versus 30 "), and are more sensitive to low mass 
cores (~ 0.5 M© versus 0.8 M© for Perseus), but they 
do not cover the clouds as completely. The sensitivity 
to density is similar, but they may be probing objects 
of somewhat lower mean density. Other differences exist 
in detail in the methods of selecting sources, so compar- 
ison to the Bolocam results will give a good idea of the 
current systematic uncertainties. iJprgensen et alj |2008) 
find more starless cores in Oph iuchus and fewer in Perseus 
than does lEnoch et al.l (|200 8b) . Calculating the lifetimes 
in the same way, we find a prestellar core lifetime of 0.55 
Myr for Perseus and 0.75 Myr for Ophiuchus. Taking the 
two clouds together yields a value of 0.61 M yr, somewhat 
larger t han the value of 0.46 Myr, obtained bv lEnoch et al.l 
(|2008aft . 

Considering the two studies, a reasonable estimate for 
the lifetime of prestellar cores above a mean density of 
about 2 x 10 4 cm" 3 would be about 0.5 ± 0.1(±0.3) Myr, 
where the uncertainty in parentheses includes the range 
over different clouds. These results are in general agree- 
ment with a third analysis of Perseus by lHatchell et alj 
(12007ft . 

In comparing our results to previous r esults , summa- 
rized in Fig. 2 of I Ward-Thompson et alj (|2007ft . our life- 
time would appear to be lower than what they found if 
plotted at our minimum density of n = 2 x 10 4 cm -3 . 
Ho wever, most cores are denser than the minimum value, 
and lEnoch et al.l (|2008bft have done a more careful anal- 
ysis, dividing the prestellar cores into two bins of mean 
density greater than or less than n — 2 x 10 5 cm~ 3 . The 
inferred lifetimes of 0.32 Myr and 0.11 Myr are quite con- 
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sistent with values found b y iKirk et alj (12005) and with 
the overall relation found bv lWard-Thompson et al.l (|2007t ) 
indicating lifetimes decreasing toward, but not yet reach- 
ing, the free-fall time (iff) as the density increases, a s 
originally suggested bv lJessop fc Ward-Thompson! (|2000T ). 
At this point, we should recall the assumptions required 
for these estimates. In particular, we assum e no depen- 
dence on core mass in deriving a lifetime. IClark et alj 
(|2007T ) sugges ted that lower mass co res collapse more 
quickly, while lHatchell k, Fuller! (|2008l ) note that, to the 



known edge -on disks is needed to evaluate this suggestion. 
In addition. ICrapsi et al.l (|2008f ) argues that half of the ob- 



contrary, a higher fraction of the more massive cores in 
Perseus contain embedded objects, suggesting faster evo- 
lution for more massive cores. The core mass distribu- 
tion for Bolocam cores with embedded sources is broader 
and flatter than th e prestellar core distribution as well 
(En och et al.ll2008bf ). also suggesting faster evolution for 
the more massive cores, though other explanations are pos- 
sible. The issues connecting the core mass function and 
the initial ma ss function of stars are explored in a num- 
ber of papers (lEnoch et al.ll2008bl lHatchell fc Fuller| [2008. 
IClark et al.ll2007l iGoodwin et al.ll2008l . lAndre et al.ll2007l 
and lSwift k Williamsll2008T) ~ 

Note that we have referred to these starless cores as 
prestellar cores because almost all with suitable molec- 
ular line data were found to be gravitationally bound 
(|Enoch et alj l2008bf ) and hence worthy of the "pr estel- 
lar" moniker, as defined by di Francesco et all ([2007) and 
IWard-Thompson et all (|2007h . 

W ith the large sample of embedded objects now avail- 
able, [Enoch~eFlil] (|2008af ) have subdivided the Class and 
Class I SEDs into early and late versions. While the SED 
evolution appears to be continuous, suitable choices of Tboi 
can produce distinctive SEDs in the mean. By counting 
the number in each bin, they found a variation in T^i 
that increases with time: T bo i = 25K + C(t/10 5 j/r) 1 - 8 , with 
C = 51 K, though the evolution must slow down to match 
the Class II lifetime. 

7.2. Mid-Life Crisis: Class I, Flat, and Class II 

The boundary between Class I and II is also problem- 
atic. As discussed in £|5.1i we will use the terminology of 
"stage" to describe the physical configuration and "class" 
to describe the SED. Mo dels with 2-D radiative tr ansfer 
by iRobitaille et alj ([2006D and ICrapsi et all (|2008ft show 
that Stage II objects with large inclination angles can have 
Class I SEDs. Stage I and II objects occupy largely dis- 
tinct regions of the various color-color diagrams but there 
is also significant over lap in some parts of these diagrams 
(Robitail le et al.| [2007). Some Stage I objects viewed rel- 
atively face-on will be classified as Class II objects using 
popular color-color diagrams such as Figure 1101 Longer 
wavelengths are helpful in distinguishing these situations; 
in particular, extended millimeter continuum emission or 
tracers of dense gas, peaked on the source in either case, 
unambiguously indicate the pr esence of an envelop e. 

Based on his grid of models, ICrapsi et al.l (|2008h argues 
that 34% of Stage II (star+disk) sources are misclassified 
as earlier types because of inclination effects. Potentially, 
all Flat sp ectrum sources cou ld be nearly edge-on Stage 
II objects (|Crapsi et al.| [2008). If correct, the effect on 
lifetimes of physical stages could be substantial; further 
comparison of the properties of Flat spectrum sources to 



jects with Class I SEDs could be part of the edge-on Stage 
II population. If correct, the lifetime for envelope dissipa- 
tion would be less than half that of t he combined l i fetime 
for Class and Class I derived above. ICrapsi et alj (|2008tl 
considered a number of methods for determining the SED 
class, including a, T^oi, and the absorption or emission in 
the silicate feature; all were subject to confusion by incli- 
nation effects. 

The most reliable indicator of how much envelope re- 
mains is the ratio of single dish to interferometer flux den- 
sity at millimeter wav elengths; these pr obe the envelope 
and disk, respectively (jCrapsi et al.ll2008h . When such ob- 
servations become available for the large number of Class 

I sources in our sample, Stage I lifetimes can be refined. 
For the three clouds with complete Bolocam data, we can 
already determine what fraction of objects we classify as 
Class I based on a can be associated with millimeter con- 
tinuum emission. Overall, 65% can be associated, larger 
than predicted by ICrapsi et all (|2008f) . If we assign all 
Flat SED class and 33% of the Class I objects to Stage 
II, the lifetime of Stage I would be 0.28 Myr. Broken 
down by cloud, the percentage that lie within one FWHM 
of a dense core is 69% for Perseus, 75% for Serpens, and 
49% for Ophi uchus. Again, these numbers differ slightly 
from those in lEnoch et alj (|2008bf ). owing to slight differ- 
ences in Class I source definition. A random distribution 
of infrared sources would produce about 5 more-evolved 
objects projected by chance on a millimeter co re for all 
three clouds taken together ([Enoch et al.ll20 08b). A sim- 
ilar comparison for the SCUBA data finds that 66% and 
53% of Class I sources in Perseus and Ophiuchus lie within 
30" of the center of a SCUBA core, and one false identifica- 
tion is predicted amon g the 66 SCUBA cores in Ophiuchus 
(porgensen et al.ll2008D . 

The large fraction of the Class I sources without 1-mm 
emission in Ophiuchus may result from the generally high 
extinction in front of the YSOs in that cloud. As dis- 
cussed earlier, high extinction will in crease the value of 
a. Th is problem has been explored by Ivan Kempen et alj 
( 2008), who mapped objects in Ophiuchus that had been 
previously classified as Class I using HCO + J = 4 — > 3 and 
C 18 OJ = 3^2 emission. They found that half (11 of 22) 
of the objects assigned to Class I by their c2d value of a 
had no HCO + emission peaked on the source. They found 
a similar result (6 of 13) for the Flat SED sources. They 
reclassified those objects as disks. The very high extinc- 
tion toward many sources in Ophiuchus had increased the 
value of a derived from observations, by amounts as high 
as 1.1. As can be seen from Table [5j extinction corrections 
decrease the number of Class I sources in Ophiuchus from 
35 to 27, using a. There is a smaller decrease for Flat 
sources, from 47 to 44. The core of Ophiuchus has even 
higher extinction than our assumed value of Ay — 9.76 
mag. Calculating separately the extinction toward Class 

II sources inside the contour of Ay = 15 in the extinction 
maps, the mean Ay is 15.7 mag inside the contour and 5.7 
mag outside. Using extinction corrections based on those 
numbers for Class I and Flat SED sources inside or outside 
that contour, the number of Class I sources is unchanged, 
but the number of Flat SED sources drops to 28. Thus, 
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very high extinction may substantially inflate the number 
of Flat (by 68%) SED and Class I (by 30%) sources as 
classified by observed a. Ophiuchus seems to be excep- 
tional in this regard, as can be seen in Figure [T2l There 
are many Flat SED sources in regions of very high extinc- 
tion, as measured by the cloud extinction maps. These 
would deredden to Class II sources. 

The objects classified as Class I by a but lacking dense 
gas tracers tended to have characteristic SEDs which 
peaked at 24 /iin and fell toward longer wavelengths. 
Similarly Class I sources without m illimeter continuum 
emission tended to h ave such SEDs (jEnoch et al.ll2008al 
iDunham et al.ll2008ft . This characteristic SED may be a 
way to identify suspicious Class I sources. Another ap- 
proach to sorting out these bord erline sources may be stud - 
ies of veiling of absorption lines (jCasali fe Matthewsf T992V 
a few sources in Ophiuchus without evidence of dense gas 
do have veilings that are genera lly associated with Class I 
objects ()Doppmann et al.ll2005l ). 



7.3. Senior Citizens: A 2D Disk Classification System 

Finally, the simple division between Class II and Class 
III sources may be obsolete . Note first that th ere is no 
minimum value for a in the[G reene et all ([1994ft scheme, 
as stars with no infrared excess, but other evidence that 
they are pre-main-sequence, can be included in Class III. 
Values as low as a = —3 are possible for the Rayleigh- 
Jeans tail of a hot star. For most of our clouds, we lack 
a complete sample of such objects, because we require a 
significant infrared excess to identify a source as a YSO, 
thus effectively setting a lower limit on a. Thus our counts 
of Class III sources, on this definition, are incomplete. 

A second point is that the classical idea that disk evo- 
lution leads to a more negative value of a is based on all 
the excess fluxes decreasing roughly together. In fact, we 
see a wide variety of SED shapes around the Class II /III 
boundary, leading us to suggest a new, two-dimensional 
classification system. 

As discussed in Cieza et"aT1 (|2007l ) and in more detail 
bv lMerm et all (|2008bft . the diversity of SEDs for star and 
disk systems can be captured better in a two-dimensional 
system. Many stars show no excess below a certain wave- 
length. We define Xtum-of / as the last wavelength with no 
detectable excess over a photosphere, as fitted to the pho- 
tometry at shorter wavelengths. Then we define a excess 
as the slope of the SED beyond Xtum-of /• Thus a excess 
characterizes the properties of the excess emission, rather 
than mixing the st ellar and dis k prope rties. 

As s hown by ICieza et afl (120071 ) and iMerin et ail 
(2008bj), classical T Tauri stars have excesses that start 
short (Xtum-of f < 2 /^m) and a tight distribution of 
otexcess, around —1. In contrast weak-line T Tauri stars 
with excesses usually have Xtum-of / > 2 /im, with a wider 
range of a excess . The spread in a e2cess values grows with 
increasing Xtum-of / until it covers the range from —3 to 
1, sim ilar to the range seen for debris disks (jChen et alj 
12005ft . 

The increase in Xtum-of f values may be caused by sev- 
eral things, including unresolved binaries, either stellar or 
substellar in mass, and inside-out photo-evaporation. The 
c2d sample includes a substantial number of "cold" disks, 
with large Xtum-of f and positive a excess , but they remain 



a small fraction of the sample, suggesting that the dura- 
tio n of this phase is sh ort. These issues will be explored 
bv IMerin et all (|2008bft . 

8. CLUSTERING 

Inspection of Figure [2] shows clearly that the YSOs are 
not spread uniformly over the cloud. The earlier classes are 
particularly tightly clustered in regions of high extinction, 
as is even more apparent in the figures in the pap ers on 
individual clouds. Furthermore, lEnoch et alj (|2007l ) have 
shown that the dense cores are also restricted to regions of 
high extinction and strongly clustered. We provide some 
simple quantification of these statements in t his section. 
We dr a w on previous studi es of clustering by ICarpenter I 
(|2000ft . lLada fe Ladal (120031) and I Allen etafl (|2007ft . 

We have categorized the environment of the YSOs us- 
ing simple criteria of volume density and number of YSOs. 
We distinguish loose (density above 1 M Q pc -3 ) and tight 
(density above 25 M Q pc~ 3 ) groups (5 < N < 35) and 
clusters (N > 35) . The number of 35 is taken from 
lLada fe Ladal (|2003ft . If the density falls below the loose 
criterion, YSOs are classified as distributed. The crite- 
rion for loose cluste ring is also taken from the review of 
lLada fe Ladal (|2003l ) as the minimum density for a cluster 
to av oid disruption by passing interstellar clouds (|Spitzerl 
119581 )- Since this criterion results in most sources within 
our clouds belonging to a loose group or cluster, we ex- 
perimented with higher density contours to find one that 
captured our intuitive notion of a tight group or cluster. 
More sophisticated analysis of the clustering behavior of 
YSOs is beyond the scope of this paper. 

The densities are computed from the surface density of 
YSOs using a nearest neighbor method, assuming local 
spherical symmetry to convert to volume density and as- 
suming a mean stellar mass of 0.5 M . All the YSOs, 
color-coded by SED class, are shown on gray-scale im- 
ages of the extinction in Figure [H The clustering prop- 
erties are discussed in mo re detail in other pape rs [e.g., 
P erseus and Ophiuchu s in lJorgensen et all (120081) : L upus 
in IMerin et all (|2008aD ; Cha II in lAlcala et alj (|200H )]. 

Here we aggregate all the clouds to discuss the total 
numbers in these different environments (Table [7J . Note 
that a source may appear more than once. For example, 
a source may lie within a tight group that itself is part of 
a loose cluster. Only 9% of the YSOs in our five clouds 
are truly distributed, while 91% are found in loose clus- 
ters and 54% are in tight clusters. Tight groups account 
for 13% and loose groups for 7%, illustrating the fact that 
most tight groups are themselves part of loose clusters. 
The low fraction of distributed YSOs is consisten t with 
previous studies (jLada et al.lll99ll . ICarpenter 1 12000ft . with 
the caveat that an older distrib uted population without 
infrared excess could be present ([Carpenter Il2000l ). 

The age of the population, measured by the ratio of 
Class I and Flat SED sources to Class II and III (Table 
[7]), depends on the environment as well. The highest ra- 
tios, hence youngest populations, are in groups, while the 
lowest ratios are found in the distributed populations, with 
clusters intermediate. The statistics for groups are heav- 
ily influenced by Perseus, which has a number of groups 
dominated by early SED classes. For clusters, the young 
cluster around NGC1333 is balanced by the older cluster 



17 



around IC348. The preponderance of more evolved sources 
in the distributed population is consistent with the idea 
that nearly all stars form in clustered environments and 
then disperse. 

The cloud crossing times, given in Table [TJ range from 
about 3 to 8 Myr. At the high end, these exceed the aver- 
age 2 Myr time over which we expect to be at least 50% 
complete to YSOs. However, Figure [2] clearly shows that 
the di stribution over P e rseus and Ophiucus is far from uni- 
form. J0rge nsen et al.l (I2008D identifies two loose clusters 



in Perseus, labeled e-NGC1333 and e-IC348 in their Ta- 
ble 6. The crossing times for these loose clusters are 2.4 
and 3.4 Myr, respectively, using the same mean speed as 
used for Table [T] Only one loose cluster can be identi- 
fied in Ophiuchus (c-Oph) and it would have a crossing 
time of 7.2 Myr. As might be expected, the YSOs fill 
this region very uneve nly (see Figure [5] and Figure 13 in 
lJ0rgensen et al.l (|2008l ) for a figure showing the boundaries 
of the cluster). For the L1688 cluster itself, the crossing 
time is 2.2 Myr. Overall, the distributions of sources are 
broadly consistent with a picture in which most stars form 
in clusters or groups and disperse at a rate reflecting the 
cloud turbulence. More detailed analysis will be needed 
before more precise statements can be made. 

The definition of loose clusters also allows computa- 
tion of a star formation efficiency for regions intermedi- 
ate between the cloud as a whole and the dense cores 
(|j0rgensen et al.|[2008D . With the revised cloud masses, 
we infer efficiencies of 9.1% and 7.7% for the e-NGC1333 
and e-IC348 loose clusters in Perseus, substantially higher 
than the 3.8% we find for the cloud as a whole. For Cha II, 
the loose cluster is nearly as large as the cloud, so the star 
formation efficienc y only increases to 3.9% for the YSOs 
(jAlcala et al.ll2008[ ). 

The highly clustered nature of star formation in these 
clouds may have implications for discussions of the Sun's 
birthplace. Both the abrupt truncation of the Kuiper Belt 
near 50 AU and the orbits of detached Kuiper Belt objects 
such as Sedna are likely the result of stellar encounters 
within the first few million years of solar system history 
(jBrown et all 12004 Kenyon & Bromley 2004). The pres- 
ence of short-lived nuclei such as 60 Fe, 41 Ca, 53 Mn, and 
26 Al in planetesimals must, in part, arise from the timely 
nucleosynt hetic input from nearby AGB stars and/or su- 
pernovae (|Wasserburg et all [20061 . The abundance and 
temporal duration of the injection of 60 Fe in particular 
has been used to argue for the formation of the Sun in a 
massive molecular cloud and dense stellar cluster complete 
with the fo rerunners of Wolf-Rayet stars and Type lb/c 
supernovae (|Loonev et al.ll2006LlBiz"zarro et al.l2 007). The 
iron yield is a strong function of metallicity, however, with 
rather lower mass stars capable of produ cing the amount 
of 60 Fe required by the meteorite data (jTachibana et alj 
2006:). Thus, the highly clustered nature of star formation 
revealed by the c2d survey and the nearby presence of early 
B stars to the Ophiuchus and Perseus clouds argues that 
such locations should also be considered as possible sites 
of the formation of the solar system. 



9. COMPARISON TO MODELS OF STAR FORMATION 
9.1. Does Turbulence Matter? 



The five clouds in this study span a small range in level 
of turbulence, as crudely probed by the mean linewidth 
of the 13 CO J = 1 — > line. We have searched for corre- 
lations between linewidth and cloud mass, star formation 
rate, depletion time, and the lifetime for Class I sources. 
None are apparent, probably because the range of turbu- 
le nce over our cloud sa mple is small. 

ISchmeia et alj (|2008t ) have used an earlier c2d delivery 
to count and classify YSOs, and they compared the na- 
ture of the clusters in Perseus, Serpens, and Ophiuchus. 
They concluded that the most centrally condensed cluster 
(Serpens) was associated with the lowest turbulence, as 
measured by the NH3 emission lines from the region of the 
cluster. In contrast, Table[T]indicates that Serpens has the 
largest linewidth averaged over the cloud, as measured in 
13 CO J — 1 — > 0. Clearly, the characterization of the tur- 
bulence can change conclusions. In addition, the earlier 
delivery was considerably more contaminated with back- 
ground galaxies. Inclusion of background galaxies would 
tend to make earlier classes look less tightly clustered. 

Some pictures of star formation invoke a very dynamic 
pictu re with protostars o r dense cores competing for gas 
(e.g.. lBonnell et alJl200l . We can assess the relevance of 
such pictures to the regions we study by computing the 
timescale for collisions between dense cores from t co u = 
{n rnre ^ r 2 v For the Serpens clu s ter A (|Harvev et al.l 
2007b), the data from [Enoch et alj (|2007l ) on n core and 
r, combi ned with a velocity disp ersion among cores of 1 
km s" 1 ([Williams fe MversfeoOOl ) yields t coU =3.1 Myr, 6 
times the Class I lifetime and 3 times the combined Class 
I and prestellar core lifetimes spent above a density of 
about 2 x 10 4 cm" 3 . For this cluster, we conclude that 
collisions of dense cores at the stage that we observe them 
are not dominant. Rather the cores we obse rve are more 
likely to evolve individually to form stars. lAndre et al.l 
( 2007) reached a similar conclusion for the main cluster in 
Ophiuchus, and even for very compact subclu sters within 
the main cluster. Also, Wal sh et al.l (2004) found very 
small (< 0.1 km s _1 ) motions of dense cores relative to 
surrounding less dense gas in a larger survey of dense cores. 
We caution, however, that the answer may be very differ- 
ent in regions forming still richer clusters, such as Orion. 

9.2. Luminosity Evolution 

We can make some comparisons of the lifetimes derived 
from the relative numbers of objects in different classes 
or stages in ^5.3 1 to mo dels of star formation. As noted 
bv lEnoch et a l. (2008b]), the timescale for prestellar cores 
to begin forming stars, once their mean densities are suf- 
ficient for detection by Bolocam surveys, is a few times 
the free-fall time at those densities, consistent w i th ear - 
lier work summarized by [Ward-Thompson et alj (|2007[ ). 
Such lifetimes are broadly consist ent with the evolution 
of ma gnetically cont rolled models ([Tassis fc Mouschoviasl 
[2005L [Adams fc Shul 120071) once the density exceeds the 
threshold for detection by B olocam. The figure from 
IWard-Thompson et al.l (|2007f ) that summarizes lifetimes 
as a function of mean density suggests that they get 
shorter, relative to a free fall time, as the density increases, 
also broadly consistent with magnetic models. 

The lifetime for the Class stage averaged over all 
clouds is about 0.16 Myr, but only 0.043 Myr in Ophi- 
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uchus. The very short lifetime often quoted for the Class 
phase (0.01 Myr) was originally derived from observa- 
tions in Ophiuchus ([Andre fe Montmerlelll994l) . but this 
now seems anomalous. If star formation in Ophiuchus has 
declined from rates that were higher about 0.1 Myr ago, 
the number of Class sources would be much lower com- 
pared to Class I sources. Such a scenario is broadly con- 
sistent with the fact that Ophiuchus also has the highest 
ratio of mass in YSOs to mass in dense gas of the three 
clouds with Bolocam data (Table 2]). Another factor that 
may have contributed somewhat to the short Class life- 
time in Ophiuchus is an overestimate of the number of 
Class I objects ( iJ7.2p because of very high extinction in 
the Ophiuchus core. 

The longer lifetime for the Class phase in most clouds 
is more consiste nt with the luminosit y distribution shown 
in Figure [T51 and [Enoch et al.l (|2008a| ). In particular, there 
is no clear evidence in our data for an early phase of very 
rapid mass accretion, which would produce luminosities 
higher than the mean values we observe. In fact, the lu- 
minosities seen in Cla ss sources are mostly similar to or 
less than predictions ([Young fe Eva ns 2005) of the evolu- 
tion in Lboi-Tboi sp ace based on simple inside-out collapse 
models (|Shulll977h . with constant accretion rates and unit 
efficiency. Mo dels with m ean accretion rates onto the star 
faster than the[Shu| (|l977f) rate are not supported by these 
observations. 

For Class I sources, m ost of the observed L^o i fall well 
below the predictions by I Young fc Evansl (|2005l ). suggest- 
ing that accretion occurs at a rate even lower than that of 
the Shu picture, which is the smallest among star forma- 
tion models. The mass in the envelope, on th e other hand, 
decrea ses about as predicted by this model (|Enoch et al.l 
2008a) . This result suggests that infall occurs at about the 
rate expected in the Shu picture, but that accretion onto 
the star occurs at a much lower rate most of the time. The 
lower rate may be partially explained if only a fraction e 
of the available mass in a dense core winds up on the star. 
As discussed earlier, comparison of the core mass function 
and the IMF suggest e « 0.3. This helps explain the low 
luminosities, but it is insufficient to explain the very low 
values that we see. In addition, the spread in L\, i of two 
to three orders of magnitude at a given T^oi range prac- 
tically demand a picture in which accretion onto the star 
is not steady. While no Class sources have L^oi < 0.1 
L Q in Figure [T3l this can be misleading; heating of a sub- 
stantial envelope just by the interstellar radia tion field can 
produ ce an apparent luminosity at this level (| Evans et al.l 
1200 ll) . The effect decreases as the envelope dissipates, but 
is still present during the Class I phase. 

The distribution of luminosities for all Class or I 
sources with millimeter emission is shown in Figure 1141 
The panels on the bottom cover 1 to 20 L Q ; there are 
three sources more luminous, with the most extreme at 70 
Lq. However the vast majority of sources have L^oi < 3 
Lq. As shown in the top panels, the trend of increasing 
numbers at lower luminosity continues down to 0.1 Lq, 
where selection effects com e into play. A more thorough 
study (|Dunham et al.|[2008t ) of the low luminosity end that 
dealt with these effects fo und 15 objects with in ternal lu- 
minosity, Li n t < 0.1 L Q . Dun ham et al . (2008) was able 
to rule out a continued rise in the luminosity distribu- 



tion to lower luminosities, which could have called into 
question whether any cores were truly lacking in internal 
sourc e s. In agreement with a n earlier analysis (|Kirk et alj 
I2007H . IDunham et ail (|2008l) found that 75% to 85% of 
cores thought to be starless before Spitzer remain so down 
to a limit of L int > 4 x 1 0~ 3 (d/140jpc) 2 L Q . As discussed 
by IDunham et alj (|2008f ) , the luminosity distribution in 
Figure [14] is quite inconsistent with constant mass accre- 
tion, which would dominate the luminosity in the early 
stages. For an object at the stellar/brown dwarf boundary 
of 0.08 Mq and a radius of 3 R©, the standard accretion 
rate from a Shu model (M± = 2 x 10 -6 M Q yr" 1 ) would 
produce L int = 1.6 L©. For the sources in Figure [T4l 59% 
have L\, i < 1.6 L Q ; Lboi generally provides an upper limit 
to Li n t for deeply embedded sources. Our results rein- 
for ce and aggravate the "luminosity proble m" first noted 
by Kenvon et alj (|1990l ). As su ggested by IKenvon et all 
( 19901) and further de v elope d by Kenvon et al. (1994) and 
IKenvon fc Hartmannl (jl995l ). non-steady accretion is the 
obvious solution, and it is nearly required by our results. 
Furthermore, it appears to begin in the very early, Class 0, 
phase. Such a scenario also invalidates any attempt to in- 
fer stellar mass from luminosity in phases where accretion 
luminosity dominates stellar luminosity. 

Non-steady accretion could be caused by a magnetic wall 
that forms in some models, resultin g in "spasmodic" accre- 
tion as the wall periodically fails ( Tassis fc Mouschoviasl 
2005). Presumably, turbulence in the infalling envelope 
could also result in variations in infall rates as well. Sim- 
ulations of star formation in turbulent clumps show vari- 
ations in ma ss accretion rates o f about one order of mag- 
nitude (e.g., lOffner et al.l l2008h . Because almost all the 
luminosity of accretion is released when material falls 
onto the star, relatively constant infall can store mass in 
a disk, producing very little luminosity until a disk in- 
stability dumps a lot of mass onto t he st ar, as in FU 
Orionis events (|Kenvon fc Hartmannl [l995). Statistical 
evidence for episodic ity extends into the Class stage 
(|Dunham et al.1 [20081 ) and there is direct evidence in one 
case for a cu rrent accretion rate m uch lower than a time 
average rate ([Dunham et al.ll2006D . 

Episodic accretion at peak rates higher than predicted 
by a Shu model would explain the few sourc es that lie 
above the tracks from lYoung fc Evansl (|2005[ ) in Figure 
[T3l Short periods of rapid accretion would presumably 
also drive strong outflows. If such episodes are more fre- 
quent in the Class period, that could explain the higher 
ratio of CO outflow lum inosity to Lboi seen in Class ob- 
jects ([Andre et al.ll2000l) . Some evidence for that picture 
might be seen in Figure [13l if one accounts for the lower 
luminosity predicted in the Class phase for constant ac- 
cretion models. 



9.3. Comparison to a Specific Model 

We can also apply our derived lifetimes to simple mod- 
els in a consistency check. Considering a Shu model ( Shu 
I1977D . the wave of infall propagates outward at the sound 
speed until it reaches the outer bound of the region that 
will collapse. This should occupy roughly half the time 
spent in the Class I phase, as the other half is occupied by 
the last shell falling inward. In this oversimplified picture, 
the feeding zone for a protostar would be given by 
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r f = Q.5at(J), (9) 

where a is the sound speed and t(I) is the lifetime for 
the Class I objects. For a purely thermal sound speed 
at T K = 10 K, 77 = 0.096t(7") pc with t{I) in Myr, or 
0.052 pc if t(I) = 0.54 Myr, or 0.042 pc using t'(I) = 
0.44 Myr . These radii are co mparable to the typical size 
of cores (|Enoch et al.ll2008bh . and the modal spacing be- 
tween YSOs in a s tudy of a large number of clusters is 
also about 0.05 pc (|Gutermuth et a l. 2008). The age we 
derive for the Class I phase, together with a Shu model, 
is consistent with the data at this level. However, core 
sizes are defined more by observational techniques than 
by evidence of a sharp boundary. In clustered regions, the 
core size is set by the overlap with nearby cores, as seen 
in projection. Sizes of cores taken from CLUMPFIND op- 
erating on the SCUBA data give somewhat smaller sizes 
(0.036 pc in Perseus and 0.018 pc in Oph iuchus), suggest- 
ing t hat resolution effects are important (jJorgensen et all 
2008). In much more crowded environments, such as the 
Orion cluster, the feeding zone could be even smaller. 

How much mass would a ccumulate in this time? For the 
inside-out collapse model (|Shulll977l ). the mass infall rate 
is given by 

M = m a 3 /G, (10) 

wit h toq = . 975 a nd a is the sound speed. As discussed 
bv lShu et al.l (|l987t) . the form of this expression is general, 
but M is constant only for certain initial density distribu- 
tions. For a temperature of 10 K, we have M = 1.6 x 10 -6 
Mq yr -1 . With an efficiency of e for matter to wind up in 
the star, M* = el. 6 x 10~ 6 Mq yr _1 , and the final stellar 
mass would be M* = 0.85e Mq i f t(I) is 0.54 Myr. For 
e = 0.3, suggested by lAlves et al.l (|2007t ). this produces a 
star of 0.25 M Q (0.20 M using t'(I) = 0.44 Myr), near 
the modal mass of the IMF. If the Flat SED objects corre- 
spond to the infall of the remnant envelope, one could add 
another 0.20 M Q in the time attributed to that class, but 
we suspect that many Flat sources do not have any sub- 
stantial envelope, as discussed above. Previous, shorter 
lifetimes for the Class I phase would have predicted lower 
final masses, so it is not trivial to understand how one 
builds any stars of mass exceeding 0.5 M Q . One could 
make more massive stars by increasing e or the sound 
speed, but the latter solution would predict even larger 
Lboi , and we already have a luminosity problem, illustrated 
in Figure[T31 Models with faster mean infall rates will have 
to explain the low mean values and broad spread of lumi- 
nosities in this figure, probably with recourse to spasmodic 
or episodic accretion. 

To test whether or not a picture of episodic mass accre- 
tion could result in a typical mass star, we consider a very 
simple model. The observed Ltoi for each source classified 
as Class or Class I according to Tt, i and associated with 
an envelope (112 sources) is turned into a mass accretion 
rate, M* by 



R is the stellar radius, assumed to be constant at 3 Rq, 
and M* is the protostellar mass, assumed to be 0.25 Mq 
for this calculation. This value of M* is chosen because 



it is half of the mean mass of 0.5 M© assumed in gl In 
reality, will grow as the protostar accumulates mass. 
However, if we assume that the average mass accretion rate 
stays constant throughout the embedded phase, the under- 
and over-estimates in M* introduced by including sources 
with masses above and below 0.25 Mq should roughly can- 
cel out. The final stellar mass accumulated, M*(t/), is 
then 

M*(t/)=5^M*i(M*), (12) 
M* 

where t(M*) is the amount of time spent in each accre- 
tion state and is given by the number of sources at each 
M+ divided by the total number of sources, multiplied by 
the Class I lifetime of 0.54 years. With a bin size in M* 
of 1.0 x 10~ 8 M yr" 1 , we calculate M*(t/) = 0.7 Mq. 
Smaller bin sizes result in identical results, while larger bin 
sizes overestimate M+(tf) since very low mass accretion 
rates are lumped together in one bin. Even though 92 of 
the 112 sources (82%) have L hol implying < 1.6 x 10~ 6 
Mq yr -1 , the end result is still a typical mass star. In 
this picture, much of the mass of stars accumulates during 
brief periods of high accretion. For this specific model and 
data sample, half the final mass is contributed during the 
0.039 Myr (7% of the Class I lifetime) represented by the 
eight highest luminosity sources. Since the duration of the 
events with the highest accretion rates is the shortest, even 
larger samples may be needed to find the rarest accretion 
e vents. 

Myers (2008) has pointed out that the dense cores do 
not have sharp boundaries, but instead inhabit "clumps" 
of lower, but still substantial, density. In this case, the final 
stellar mass may be determined by a competition between 
free fall onto the protostar and dispersal of the clump gas. 
More massive stars could be formed in this way. If the 
infalling material from the clump has sufficient opacity to 
make the SED fall into Class I, the time limit on this pro- 
cess would still be the lifetimes derived above, and those 
would limit the final stellar mass to values similar to those 
derived above. 

To compare more specifically to theoretical models, it 
is necessar y to translate the mod els into simulated ob- 
servations. lYoung k. Evans! (|2005l) did that for the stan- 
dard Shu model (Shu 19771 ). assuming unit efficiency and 
spherical radiative transfer. They predicted that the ob- 
served Class phase would be short, about 0.05 Myr, 
terminated by the collapse of the first hydrostatic core 
as the central source shrank to 3 Rq and accretion lu- 
minosity became important, and significant emission in 
the mid-infrared began to be produced. That lifetime was 
very uncertain because it depended on assumptions about 
the first hydrostatic core, so the time may be longer; al- 
ternatively, episodic accretion could move sources back 
and forth across the Class 0/1 boundary. The Class 
sources with high L\> i would need to have thicker en- 
velopes or be viewed edge-on through a non-spherical enve- 
lo pe. The transition from Class I to Class II in the models 
of lYoung fc Evansl (|2005t ) depended strongly on the mass 
of the initial envelope, ranging from 0.1 to 0.4 Myr for 
masses of 0.3 to 3 Mq. This transition occurred near, but 
not at, the end of infall. The time spent in the range of 
Tboi we might attribute to the Flat stage was very short, 
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reaching 0.1 Myr only for the model with a mass of 3 Mq. 
Clearly, the effects of non-spherical models and non-steady 
accretion will have to be considered to bring observations 
and models into better congruence. 

10. FUTURE WORK 

While the c2d study has provided answers to many ques- 
tions, much work remains to be done. It will be important 
to incorporate results from the Taurus and Gould Belt 
surveys. Spectral types for the remaining visible sources 
are necessary, both to weed out remaining contaminants 
and to estimate masses of the forming stars. Searches for 
binaries would be very valuable, especially for the cold 
disks. More complete lists of PMS stars without infrared 
excess are needed in most clouds. Near-infrared surveys 
that go much deeper than 2MASS would better match the 
sensitivity of the Spitzer data, aiding greatly in discrimi- 
nating against background galaxies. Low-noise millimeter 
continuum maps of Lupus and Cham II are needed to put 
those clouds on the same footing as the northern clouds. 
Episodic accretion could leave signatures in the envelope 
chemistry or in outflows, so studies of multiple molecular 
lines are needed for the embedded sources. 

11. SUMMARY 

The data from five large clouds studied by the Spitzer 
c2d legacy project have been combined and reanalyzed in 
a consistent way. The five clouds combined contain 1024 
objects we classify as YSOs. We estimate that there are 
no more than about 50 (5%) contaminating background 
galaxies in the sample and probably many fewer because 
of human examination. We estimate that the sample is 
complete down to L^oi ~ 0.05 L Q over a wide range of 
SED types. This limit arises from increasing contami- 
nation by background, star-forming galaxies at low flux 
densities rather than sensitivity limits. With better dis- 
crimination against galaxies, the fundamental luminosity 
limit could be lowered to about 10 -3 L©. 

The five clouds are forming about 260 M© of stars per 
Myr. Using cloud masses determined from c2d extinction 
maps, the star formation efficiency over the last 2 Myr is 
3% to 6%, depending on the cloud. The time to deplete the 
current cloud material at the current rate exceeds most es- 
timated cloud lifetimes, but final efficiencies could be 15% 
to 30% if the clouds form stars at the current rate for 10 
Myr. That would require continued creation of dense cores 
from less dense cloud material at a rate that replenishes 
the current stock in about 2 Myr. 

The star formation surface density is about 20 times 
larger than would be predicted by relations u sed for ex- 
tragalactic star formation, whethe r non-linear (jKennicuttj 
119981) or linear (jBigiel et all [20081 in gas surface density. 
The result s are more con sistent with relations found for 
dense gas (|Wu et al J 120051 ) . 

For three clouds, comparisons to dense gas, as mapped 
by continuum emission from dust, are possible. The effi- 
ciencies are much higher when measured against the avail- 
able mass of dense gas. Depletion times for the dense gas 
are about 1.8 Myr, similar to timescales for cluster forma- 
tion. However, star formation is still slow compared to a 
free-fall time, even in the dense gas. 



We use the spectral index in the near-infrared to mid- 
infrared to place YSOs into traditional SED classes. Com- 
bining the number counts with an estimate of 2 Myr for the 
lifetime of Class II YSOs, and assuming a continuous flow 
through the classes, we derive lifetimes for prestellar, Class 
0, Class I, and Flat SED phases of 0.46, 0.16, 0.54, and 
0.40 Myr, respectively. If corrections to flux densities are 
made for measured or estimated extinctions before com- 
puting evolutionary indicators, the lifetimes become 0.47, 
0.10, 0.44, and 0.35 Myr for prestellar, Class 0, Class I, 
and Flat SED classes. These lifetimes should be viewed as 
median lifetimes or half-lives, rather than a lifetime that 
applies to all objects. The prestellar lifetimes refer to cores 
whose mean density is at least 2 x 10 4 cm~ 3 and more typ- 
ically above 10 5 cm -3 . Life times of cores with lower mea n 
densities are clearly longer (|Ward-Thompson et al.|[2007t ). 
and may vary fro m cloud to cloud [cf. the Pipe nebula 
(lLada et al.ll200l ]. 

The connection between the traditional classes and ac- 
tual physical stages is uncertain, particularly in regions of 
very high extinction, such as the core of Ophiuchus. Class 
II objects behind heavy extinction can be classified as Flat 
or even Class I objects. Our statistics after extinction 
corrections illustrate the uncertainties (typically 10% to 
20%), but the effects can be la rger. As discussed in more 
detail by iMerm et al.l (|2008bl ). the traditional evolution 
from Class II to Class III does not capture the diversity 
in SEDs found in our survey. We suggest a two-parameter 
classification as an intermediary to more physical interpre- 
tations. 

Star formation is far from uniformly distributed over 
the clouds. YSOs are highly concentrated in regions of 
high extinction and the younger YSOs are even more con- 
centrated in regions where millimeter continuum emission 
indicates the presence of dense cores. The majority (54%) 
of YSOs are in tight clusters (at least 35 YSOs and a den- 
sity above 25 M© pc~ 3 ), and 91% lie within loose clusters 
(at least 35 YSOs and a density above 1 M© pc~ 3 ). Only 
9% are truly distributed. 

We find no obvious trends with the degree of turbulence 
in the clouds as a whole, but the range is very small. Esti- 
mates of collision times between dense cores suggest that 
they arc unlikely to collide during their remaining lifetime, 
suggesting that individual cores may retain their identity 
and form a small number of stars. Thus our data are con- 
sistent with a mapping of cor e mass function onto stellar 
mass function, as discussed bv lEnoch et alj (|2008b| ). How- 
ever, denser star-forming clusters, such as Orion, may be 
quite different. 

Comparison to a simple inside-out collapse model sug- 
gests that a star of mass 0.25 M© could form during the 
lifetime of the Class I stage, assuming an efficiency of 0.3 
for infalling mass to wind up in the star. However, we 
see a spread of at least three orders of magnitude in L^i 
in the embedded phases, with most YSOs below the lu- 
minosity expected in an inside-out collapse model. The 
data strongly indicate that accretion onto the forming star 
must be transient, with very large fluctuations. Thi s result 
confirms the conclusion of iKenvon "fc~ Hartmann ( 1995) 
and extends it to earlier stages. A very simple model 
of episodic mass accretion applied to our sample suggests 
that a typical star could assemble half its mass in a few 
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episodes of high accretion, corresponding to the high lumi- 
nosity sources in our sample. These episodes could occupy 
only about 7% of the full Class I lifetime. 
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Table 1 
Facts about Clouds 



Cloud 


Solid angle 


Distance 


Area 


Av 


Mass a 


(n) b 


t (cross) 


Rcfs 




(deg 2 ) 


(PC) 


(pc 2 ) 


(km s- 1 ) 


(M ) 


(cm -3 ) 


(Myr) 




Cha II 


1.038 


178 ± 18 


10.0 ±2.0 


1.2 


426 ± 86 


345 


3.7 


1, 2 


Lupus 


3.101 


150 ± 20 c 


28.4 ±6.5 


1.2 


816 ± 188 


381 


4.7 d 


3,4 


Perseus 


3.864 


250 ± 50 


73.6 ±29.4 


1.54 ±0.11 


4814 ± 1925 


196 


7.8 


5, 6 


Serpens 


0.850 


260 ± 10 


17.5 ± 1.4 


2.16 ±0.01 


2016 ± 155 


707 


2.7 


7, 6 


Ophiuchus 


6.604 


125 ±25 


31.4± 12.6 


0.94 ±0.11 


2182 ± 873° 


318 


8.4 


8, 6 


Total 


15.457 




160.9 ±51.9 




10254 ± 3228 


389 







"The masses are computed from extinction maps with 240" resolution based on c2d and 2MASS data, except for Ophiuchus, which uses the 
270" resolution map; the mass refers to the area with Ay > 2 mag, and the uncertainty reflects only the distance uncertainty. 

''The mean density of the cloud, calculated from the mass and the surface area, assuming a spherical cloud. For Lupus, the value is an average 
over the three clouds. 

c The Lupus III cloud is at 200 ± 20 pc. This is accounted for in the total area and mass. 

''This is the crossing time for Lupus III, the largest subcloud. The time for Lupus I is 3.6 Myr and for Lupus IV, it is 1.9 Myr. 
e This mass excludes Ophiuchus-North, a disconnected piece of the northern streamer. 
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Table 2 

Number of Objects of Various Types in Catalogs 



Category 


Cha II 


Lupus I 


Lupus III 


Lupus IV 


Perseus 


Serpens 


Ophiuchus 


Full Archive 


265607 


410067 


480766 


177727 


777484 


377456 


1629665 


High Reliability 


41787 


47597 


101507 


33340 


58312 


104098 


211638 


Three Bands 


47952 


67494 


113515 


47367 


70071 


91555 


207496 


Stars 


22958 


30371 


56813 


21494 


25753 


57784 


106964 


YSOc 


29 


20 


79 


12 


387 


262 


297 


GALc 


298 


337 


299 


79 


826 


208 


842 


Other 


5036 


6539 


9859 


3026 


11529 


9807 


207496 



Table 3 

Numbers, Densities, Star Formation Rates 



Cloud 


N(YSOs) 


N/fi 
(deg" 2 ) 


N/Area 
(pc" 2 ) 


SFR 

(Mo Myr- 1 ) 


SFR/Area 
(M Q Myr^pc" 2 ) 


Notes 


Cha II 


26 


25 


2.6 


6.5 


0.65 


Alcala et al. (2008) 


Lupus a 


94 


30 


3.3 


24 


0.83 


Merin et al. (2008a) 


Perseus 


385 


100 


5.2 


96 


1.3 


Lai et al. (2008) 


Serpens 


227 


267 


13.0 


57 


3.2 


Harvev et al. (2007t>) 


Ophiuchus 


292 


44 


9.3 


73 


2.3 


Allen et al. (2008b) 


Total 


1024 


66 


6.4 


256 


1.6 





a A sum over the three Lupus clouds of the values for each cloud. 



24 



Table 4 

Efficiencies and Depletion Timescales 



Cloud 


(M(cloud)+M,) 


AU/M (dense) 


t dep (cloud) 
(Myr) 


t dep (dense) 
(Myr) 


SFRff 


Notes 


Cha II 


0.030 




66 




0.028 


Alcala et al. (2008) 


Lupus a 


0.054 




35 




0.050 


Merm et al. (2008a) 


Perseus 


0.038 


0.69 


50 


2.9 


0.049 


Lai et al. (2008) 


Serpens 


0.053 


1.2 


35 


1.6 


0.036 


Harvev et al. (2007b) 


Ophiuchus 


0.063 


3.3 


30 


0.6 


0.064 


Allen et al. (2008b) 


All Clouds b 


0.048 


1.2 


40 


1.8 


0.040 





a A sum over the three Lupus clouds of the values for each cloud. 

6 For all but SFRff, this number is calculated by adding all clouds with the relevant data together; for SFRff, it is the average over all clouds 
of the individual values. 



Table 5 

Number of YSOs by Cloud and Class 



Cloud 


I 


Flat 


II 


III 


t(I) 


t(F) 


Reference 












(Myr) 


(Myr) 




Cha II 


2 


1 


19 


4 


0.21 


0.11 


Alcala et al. (2008) 


Lupus a 


5 


10 


52 


27 


0.19 


0.38 


Merm et al. (2008a) 


Perseus 


87 


42 


225 


31 


0.77 


0.37 


Lai et al. (2008) 


Serpens 


36 


23 


140 


28 


0.51 


0.33 


Harvev et al. (2007b) 


Ophiuchus 


35 


47 


176 


34 


0.40 


0.53 


Allen et al. (2008b) 


Total 


165 


123 


612 


124 


0.54 


0.40 




After E.C. b 










t'(I) 


t'(F) 


A v Used 


Cha II 


1 


2 


18 


5 


0.11 


0.22 


3.95 


Lupus a 


5 


7 


54 


28 


0.19 


0.26 


2.91 


Perseus 


76 


35 


244 


30 


0.62 


0.29 


5.92 


Serpens 


32 


22 


140 


33 


0.46 


0.31 


9.57 


Ophiuchus 


27 


44 


179 


42 


0.30 


0.49 


9.76 


Total 


141 


110 


635 


138 


0.44 


0.35 





a A sum over the three Lupus clouds of the values for each cloud. 
'After correction for extinction, as described in £13.31 



Table 6 

Dependence of Class on Definition 



Method 


Class 


Class I 


Flat 


Class II 


Class III 


Notes 


Catalog 




165 


123 


612 


124 


fit to 2-24 


IRAC 




145 


104 


580 


187 


fit to 3.6-8 


De-reddened 




141 


110 


635 


138 


fit to 2-24 


T 6o ;(obs) 


40 


156 




753 







T fooi (dered) 


24 


125 




223 


577 




T b oi{ohs) 


44 


88 










with mm core 


T bo i{dered) 


26 


95 




11 





with mm core 



"The Catalog values of a are from the c2d catalog, based on a least-squares fit to any available data from 2 to 24 /im. IRAC values are based 
only on IRAC wavelengths. De-reddened values of a or T(, ;are based on extinction estimates from spectral types where available, assumption 
of a K7 spectral type otherwise for Class II and III sour ces. For Flat and Cl ass 0/1 sources, the mean extinction toward Class II sources in 
that cloud was assumed. The last two lines are based on Enoch et all J2008a), who require a detection with Bolocam and thus refer only to 
the three clouds with such data. 
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Table 7 

Sources by Class and Environment 



Environment 


Class I 


Flat 


Class II 


Class III 


Total 


I+F/II+III 


Distributed 


11 


7 


43 


32 


93 


0.24 


Loose Group 


28 


5 


30 


7 


70 


0.89 


Tight Group 


34 


25 


64 


8 


131 


0.82 


Loose Cluster 


127 


112 


559 


131 


929 


0.35 


Tight Cluster 


90 


79 


322 


63 


554 


0.44 
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APPENDIX 
TABLES OF YSOS 

Table IA8I contains the list of YSOs analyzed in this paper. The columns include a running index number, which is 
consistent throughout these tables, the JHK flux densities from 2MASS, flux densities from the four LRAC bands, and 
flux densities from the two shortest wavelength bands of MIPS. Additional flux densities from observations with other 
teles copes, at wav e lengths ranging from 0.36 /x m to 1.3 mm are p r esented in T ables 1A91 ETOl and lAlll The d ata are taken 
fromlAlcala et all (12008TI. iBrede et all (l2008lh Ipavis et alJ (ll999h.lDiupvik et all (12009). lEnoch et al l (T2006). E noch et all 
(120071). lHatchell et all (120051). iHatdiell et all (12007D. iHenning et all (119931) Uohnstone et all (l2000f). iKirk et alj (1200"" 
iMerin et all (|2008aD. IStanke et all (|2006f). iTachihara et all (|2007l).IWu et all (|2007l). andlYoung etall (|2006h 



Table |A"T2"1 contains the derived properties for each source. The columns list, in order, the running index number, the 
name of the cloud containing the source, the official c2d source name, which includes the position information in J2000 
coordinates, the c2d source classification, and the values for a, Tboi, and Lboi, calculated from the observed flux densities. 
Following are the same three quantities (a', T b ' oi , and L' bo[ ), calculated after correcting the flux densities for extinction, as 
described in £13.31 The nex t column ind i cates whether there is ev i dence for an envelo pe, as traced by extended millimeter 
continuum emission (e.g., En och et alj I2006T lYoung et alj [20061 lEnoch et all 12007) . The last column has information 
on the reliability of the Tboi and Lboi values. "G" indicates existing photometry provides sufficient spectral coverage 
for reliable calculations; the uncertainties in both quanti ties are dominated by errors introduced by incomplete, finite 
sampling of the source SED and are typically 20 — 60% (jEnoch et al.ll2008at iDunham et all 120081 ). "L_I" indicates the 
source is associated with an envelope and thus likely an embedded source, but lacks spectral coverage between 24 /xm and 
1 mm. The calculated values of Lboi and Tboi are thus considered lower and upper limits, respectively. "LJI" indicates 
the source is not associated with an envelope and has an observed a < —0.3, thus is not likely an embedded source, but 
no extinction correction could be derived based on available data, usually because of a lack of sufficient coverage in the 
near-infrared. The calculated values of Lboi and Tboi are thus considered lower limits. 

The full tables are quite long and available only electronically. We provide samples of the tables in the print version. 




Fig. 1 — 



Images of all of the clouds on the same absolute size scale. 



The color code is blue (4.5 (im), green (8.0 ^tm), and red (24 ^m). 




Fig. 2. — Gray scale images of the extinction of all of the clouds on the same absolute size scale as shown by the 3 pc scale bar. All YSOs 
are plotted with the following colors red (I), green (flat), blue (II), and purple (III). The gray scale is linear from Ay = 1 to 25 mag. 
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Fig. 3. — A histogram of the number of objects in bins of the log of the probability that the source is a galaxy. The combined data for the 
five large clouds are represented by the solid line and the results for the SWIRE data (degraded to simulate our data) are shown with dotted 
lines. The vertical dashed line at logP(gal) = —1.47 is the threshold used to separate YSOc from Gale (galaxy candidates). 
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Fig. 4. — The s urface density of star formation is plotted versus the mass surface density. The solid line marked K98 shows the relation from 
Kcnnicutt (1998), and the filled circles show our data. The three Lupus clouds are plotted separately. The open square shows the average for 
all our clouds. The line labeled Bigiel08 shows the linear relation from iBigiel et all (120081) (solid over the range o f their study an d dotted as 
extrapolated to higher surface densities). The line labeled Wu05 shows the relation for dense gas traced by HCN (Wu ct al. 2005) (solid over 
the range studied in dense cores in our Galaxy and dashed as extrapolated to lower surface densities). 




Fig. 5. — Pie charts for each cloud and for the total, showing the percentage of sources in each SED class. Red is Class I, green is Flat, 
blue is Class II, and purple is Class III. 
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Fig. 6. — The value of a is plotted versus T bo i for each source with well-determined values. The color code is based on the Lada class of 
each source, as defined by a, with Class I plotted as red, Flat as green, Class II as blue, and Class III as purple. Filled circles indicate sources 
associated with envelopes as traced by millimeter continuum emission, while plus signs indicate sources with no such associations. The right 
panel shows a' and T^ ol calculated after corrections for extinction were applied to the observations, as described in i|3.3l The vertical dashed 
lines show the boundaries between classes, as defined by Chen ct al. (1995). 
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Fig. 7. — The distribution of Tj, ; for sources in each of the Lada classes, as defined by a, in linear T(, ; bins. The color code is the same 
as in Figure [5] (Class I plotted as red, Flat as green, Class II as blue, and Class III as purple). The right panel shows the same distributions 
after corrections for e xtinction were app lied to the observations, as described in i|3.3l The vertical dashed lines show the boundaries between 
classes, as defined bv lChen et al.l 1 119951) . 
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Fig. 8. — The c2d source type is shown versus T;, ; for each source with well-determined values. The color code is the same as in Figure [6] 
(Class I plotted as red, Flat as green, Class II as blue, and Class III as purple). Also as in Figure [6] filled circles indicate sources associated 
with envelopes as traced by millimeter continuum emission, while plus signs indicate sources with no such associations. The right panel shows 
the c2d source type versus T^ o[ , calculated after corrections for extinction were applied to the observations, as described in i|3.3l The vertical 
dashed lines show the boundaries between classes, as defined by Chen et al. (1995). 
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Fig. 9. — The c2d source type is shown versus a for each source with well-determined values. The color code is the same as in Figure [6] 
(Class I plotted as red, Flat as green, Class II as blue, and Class III as purple). Also as in Figure [6] filled circles indicate sources associated 
with envelopes as traced by millimeter continuum emission, while plus signs indicate sources with no such associations. The right panel shows 
the c2d source type versus a' , calculated after corrections for e xtinction were applie d to the observations, as described in i|3.3l The vertical 
dashed lines show the boundaries between classes, as defined by Greene ct al. (1994). 
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[5.8]-[8.0] 

F ig, 10. — The color-color diagram for the 4 IRAC bands. The diagram on the left is constructed from the observed photometry in Table 
IA8I while the one on the right is constructed after approximate corrections for extinction. The color code is the same as in Figure [6] (Class 
I plotted as red, Flat as green, Class II as blue, and Class III as purple) with the addition that Class sources are plotted in black. Also 
as in Figure [6] filled circles indicate sources associated with envelopes as traced by millimeter c ontinuum emission , while plus signs indicate 
sources with no such associations. The box indicates the area identified with Class II sources by I Allen etaLI J2004T) . 



35 




Fig. 11. — The color-color diagr am u sing 3 of the 4 IRAC bands and the first MIPS band. The diagrams on the left are constructed from 
the observed photometry in Table IA8I while the ones on the right are constructed after approximate corrections for extinction. The color 
code is the same as in Figure [TOl (Class plotted as black, Class I plotted as red, Flat as green, Class II as blue, and Class III as purple). 
Filled circles indicate sources associated with envelopes as traced by millimeter continuum emission, while plus signs indicate sources with no 
such associations. The boxes in the upper two panels indicate the areas identified with Class Ill/stellar, Class II, and Class 0/1 sources by 
IMuzerolle et al.l l|2004n . going from lower left to upper right. The botto m two panels overlay the areas filled by Stage III, II, and I sources 
moving from lower left to upper rig ht, taken from IRobitaille et al.l ll2006h . 
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Fig. 12. — The extinction for each YSO in Perseus, Ophiuchus and Serpens versus its value of a. The extinction for each source is extracted 
at its position in the map of extinction based on background st ars. The extinction vec t or (re d arrow) shows the effect on a assuming that 
half the extinction lies in front of the source. The opacity law of Wcinga rtner &: Draind 1 1200 II ) with Ry = 5.5 is used. The gray-scale shows 
the surface density of points. Ophiuchus has a distinct population of flat spectrum sources at an Ay of around 30, which could be reddened 
from the main locus of the Class II sources at l ow extinctions. The vertical dotted lines show the boundaries between Class II, Flat, and 
Class I sources as defined bv lGreene eTaLl jl994h . 
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Fig. 13. — The bolometric luminosity (Lf, i) is plotted versus the bolometric temperature (Tt, ;) for all sou rces with sufficient d ata to 
properly constrain the values. Similar plots for the embedded sources only and for the individual clouds are in lEnoch et a l. (2008a). The 
color code is the same as in Figure [6] (Class I plotted as red, Flat as green, Class II as blue, and Class III as purple). Also as in Figure[5] filled 
circles indicate sources associated with envelopes as traced by millimeter continuum emission, while plus signs indicate sources with no such 
associations. The three thick lines are model tracks for initi al core masses of 0.3, 1.0, and 3.0 Mq, collapsing according to a Shu inside-out 
collapse model in which all mass becomes part of the star (Young & Evans 2005); they end when all infall stops, but the later stages are 
highly uncertain. The three thin lines are model tracks with initial core masses of 1.0, 1.8, and 3.0 Mq, collapsing with an accretion rate tha t 
decreases exponentially with time and efficiencies such that the final stellar masses are 0.5, 0.3, and 0.5 Mq, respectively (Myers ct al. 1998). 
The right panel shows L' bol and T' hol calculated after corrections for extinction were applied to the observations, as described in £|3.3I The 
vertical dashed lines show the boundaries between classes, as defined bv lChen et al.l 1 119951) . The heavy dashed line on the left is the ZAMS 
(D'Antona & Mazzitclli 1994) from 0.1 to 2 M . 
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Fig. 14. — The number of embedded sources (associated with envelopes) in luminosity bins. The lower plots show all sources up to 20 
Lq. There are three sources with higher L boi , with L bo i ~ 70 Lq as the largest. The top panels show a blow up of the range below 1 Lq. 
The panels on the right show the histograms obtained after correcting for extinction. The vertical dashed line indicates the level of L bo i that 
can easily result just from heating from the interstellar radiation field. This value can vary by factors a few, depending on the mass of the 
envelope and t he strength of the rad iation field. There is a selection bias against lower luminosities. This low luminosity regime is studied in 
more detail by [Dunham et all (t2008l 'l. 
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Table A8 

2MASS and Spitzer Flux Densities, in milli-Jy, of YSOs in the c2d Clouds 
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51 ± 1.5 


78 ± 2.4 


86 ± 2.5 


73 ± 4.0 


67 ± 3.3 


61 ± 2.9 


57 ± 2.7 


51 ± 4.7 


40 ± 5.7 


6 


12 ± 0.35 


1.0 ± 0.29 


8.7 ± 0.25 


5.9 ± 0.29 


4.9 ± 0.24 


3.9 ± 0.20 


3.1 ± 0.15 


0.82 ± 0.21 




7 


160 ± 4.4 


270 ± 7.5 


300 ± 8.2 


290 ± 15 


250 ± 13 


260 ± 13 


330 ± 16 


380 ± 36 


170 ± 19 


8 


1.5 ± 0.060 


4.2 ± 0.12 


6.3 ± 0.21 


7.1 ± 0.35 


6.8 ± 0.34 


6.7 ± 0.33 


7.0 ± 0.33 


8.1 ± 0.76 




9 




2.7 ± 0.12 


38 ± 1.1 


210 ± 10 


440 ± 26 


640 ± 31 


700 ± 34 


3600 ± 340 


6200 ± 600 


10 


1200 ± 29 


1800 ± 55 


1700 ± 46 


800 ± 41 


450 ± 24 


390 ± 19 


320 ± 15 


260 ± 24 





Table A9 

0.36 — 0.96 fiM Flux Densities, in milli-Jy, of YSOs in the c2d Clouds 
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Table A 10 

3.4 - 100 fiM Flux Densities, in milli-Jy, of YSOs in the c2d Clouds 
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Table All 

160 - 1300 fiM Flux Densities, in milli-Jy, of YSOs in the c2d Clouds 
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Table A12 
Properties of YSOs in the c2d Clouds 







Spitzer 










Observed 




Extinction Corrected 








Source Name 


c2d 






Tbol 


Lbol 




Tbol 


Lboi 






Index 


Cloud 


(SSTc2d +) 


Classification 


A v a 


a 


(K) 


(L©) 


a 1 


(K) 






Quality 


1 


Cham II 


J124825. 74-770636. 5 


YSOc_star+dust (MP 1 ) 


4.4 


-2.19 


2200 


11 


-2.59 


3000 


28 


N 




2 


Cham II 


J125230.66-771513.0 


YSOc_star+dust(MPl) 


4.3 


-2.38 


2100 


1.4 


-2.71 


2700 


3.2 


N 




3 


Cham II 


J125317.23-770710.7 


rising 


10.5 


-0.72 


550 


38 


-0.14 


1400 


72 


N 




4 


Cham II 


J125342.86-771511.5 


YSOc_red 


4.0 


0.65 


120 


0.49 


0.21 


190 


0.59 


Y 




5 


Cham II 


J125633.66-764545.3 


YSOc_star+dust(IR2) 


2.4 


-1.22 


2200 


0.21 


-1.35 


2900 


0.36 


N 




6 


Cham II 


J125658.68-764706.6 


YSOc_star+dust(IR2) 


0.0 


-1.84 


>4500 


>0.067 


-1.84 


>4500 


>0.067 


N 


L. 


7 


Cham II 


J125711. 77-764011. 3 


YSOc_star+dust(IRl) 


3.3 


-0.85 


1800 


0.86 


-1.06 


3100 


1.8 


N 




8 


Cham II 


J125806.78-770909.4 


YSOc_star+dust(IRl) 


5.0 


-0.93 


1400 


0.013 


-1.09 


1900 


0.020 


N 




9 


Cham II 


J125906.58-770739.9 


YSOc_red 


4.0 


0.68 


220 


1.9 


0.49 


270 


2.4 


Y 




10 


Cham II 


J125926.45-774708.4 


YSOc_star+dust(IR4) 


3.9 


-1.78 


2400 


3.7 


-2.01 


3300 


9.8 


N 





"Value of Av used for dereddening, as explained in the text. 

^Indicates whether or not the source is associated with an envelope as traced by extended millimeter emission (based primarily on the sample 

of Enoch et al. [2008]). 
indicates the quality of the calculated Lj, ; and T(, D j values. 

"G" : existing photometry provides sufficient spectral coverage for reliable calculations. Other flags mark sources without sufficient coverage: 

"L_I": the source is associated with an envelope and thus likely an embedded source, but lacks spectral coverage between 24 /im and 1 mm. 

The calculated values of L 00 [ and T(, ; are thus considered lower and upper limits, respectively. 

"L_II": the source is not associated with an envelope and has an observed a < —0.3, thus is not likely an embedded source, but no extinction 
correction could be derived based on available data, usually because of a lack of sufficient coverage in the near-infrared. The calculated values 
°f ^bol an d Tbol arc thus considered lower limits. 

"O": Other. These three sources were added to the list of YSOs based on the searches for embedded protostars presented by j0rgensen et al. 
(2007) (for Perseus) and j0rgcnscn et al. (2008) (for Ophiuchus), but they are not associated with envelopes according to Enoch et al. (2008). 
They are only detected at one wavelength (24 fim), thus calculations of a, T(, ;, and L\y i arc not possible. They are counted as Class I objects 
when classifying by a and a'; they are not counted at all when classifying by Tf, i and T\y i . 



